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FOREWORD 

This report is one of a series of technical documents prepared for the 
Grand River Basin Water Management Study. The project described herein was 
undertaken by the Grand River Hydrology Subcommittee at the request of the 
Grand River Implementation Committee. 

The material contained in this report is primarily technical support 
information and in itself does not necessarily constitute policy or 
management practices. Interpretation and evaluation of the data and 
findings in most cases, cannot be based solely on this report, but 
should be analysed in light of other reports produced within the 
comprehensive framework of the overall study. Questions with respect to 
the contents of this report should be directed to the Co-ordinator of te 
Grand River Water Management Study, Water Resources Branch, Ministry of the 
Environment, 135 St. Clair Avenue West, Toronto. 
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1.0 BACKGROUND 



Technical Report No. 35 is one of a series of five technical 
reports dealing with flood flows and flood damages in the Grand River 
basin. These reports are: 

Technical Report No. 35 Flow Frequency Analysis 

Design Floods 



No. 


35 


No. 


36 


No. 


37 


No. 


38 


No. 


39 



Application Of The HEC-5 Model 
Reservoir Operations 
Flood Damages 

The overall objective of the flood flows study is to produce 
stage-flow-frequency-damage relationships at all flood prone urban 
centres: 

a) under natural conditions; 

b) under regulation by the existing reservoir system, and as 
affected by existing dyke systems; 

c) under regulation by proposed reservoirs, and as affected by 
proposed dyke and channel improvement schemes. 

In order to satisfy this objective, several steps were undertaken: 

In Technical Report No. 35 (this report): 

1) the natural flood hydrograph was estimated for maximum annual 
spring freshet events at long-term gauge stations where flows 
are regulated by the existing reservoir system; 

2) the flow-frequency relationship was established at: 

a) the long-term stations using the natural peak flows 
established in step 1. 

b) at long-term stations where flows have not been 
regulated. 

In Technical Report No. 36: 

3) the hydro! ogical processes of the watershed were modelled; 
the 5, 10, 25, 50 and 100-year flood flows at all flood prone 
urban centres were obtained by simulating flood events whose 
peaks at stations with long-term record matched the natural 
flow-frequency relationship established in step 2. 

In Technical Report No. 38: 

4) the operation of the major reservoirs during high flow events 
was modelled; the effects of the reservoir systems, existing 
and proposed, on the 5, 10, 25, 50, and 100-year and the 
regional storm flood events were determined. 

In Technical Report No. 39: 

5) the frequency-damage relationships were estimated at each 
flood prone urban area for existing conditions and with 
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proposed dyking, channel improvement, and flood proofing 
schemes in place. From these relationships, average annual 
damages were calculated. 

While the study was undertaken with the overall objective in mind, the 
results of the flood frequency analysis have potentially wide 
application and a separate reporting of the results is justified. 

The purpose of this report is to present the results of the flood 
frequency analysis for gauge stations of long term record within the 
grand river basin. 



2.0 Sources Of Data 



2.1 Streamflow Data 



Several streamflow gauges have been in operation in the Grand River 
basin by the Grand River Conservation Authority and/or Water Survey 
Canada for at least thirty years. Information on these gauges is 
summarized in Table 2-1. The locations of the gauges are shown in 
Figure 2-1. 

Table 2-1 - Long Term Streamflow Gauges In The Grand River Basin 



of 



WSC Number 



Location 



Period Of 
Daily 



Record 
Continuous 



Reg. Of Flow 



02GA016 Grand River 

below Shand Dam 

02GA003 Grand River at 
Gait 



1913-29 1930 



Shand dam after 
1950- 1942 

Shand dam after 
1942; Conestogo 
Dam after 1957; 
Guelph Dam after 
1975 



02GB001 


Grand River at 
Brantford 




1948- 


Shand Dam after 
1942; conestogo 
dam after 1957; 
Guelph Dam after 
1975. 


02GA015 


Speed River 
below Guelph 


1950-60 


1961- 


Guelph Dam after 
1975 


02GA018 


With River at 
New Hamburg 


1950-64 


1965- 




02GA010 


Nith River at 
Canning 




1947- 





High flows in the Grand River basin are normally generated by one of 
three types of storms: spring rain on snow, summer thunderstorms, or 
late summer topical storms. The most frequent floods are caused by 
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spring rain on snow. Thunderstorms are not frequently wide enough in 
areal extent to cause an annual maximum flood on the Grand River or 
its major tributaries. Tropical storms do not frequently have a major 
effect on Southern Ontario. In order to carry out frequency analysis 
on a homogeneous sample, the annual spring maximum flows (those 
occurring in February through May) were first considerd. This sample 
precludes the use of the flood event caused by Hurricane Hazel. The 
effect of including Hurricane Hazel in the sample was also 
investigated (see Appendix F). 

2.2 Effect Of Reservoirs 

Three major reservoirs have been constructed in the last forty years 
for the purposes of flood control and summer flow augmentation: Shand 
Dam and Bel wood Lake (1942), Conestogo Dam and Lake (1957), and Guelph 
Dam and Lake (1975). These reservoirs, having a maximum flood storage 
(in February) of 120,000 acre-feet, have a significant controlling 
effect on the flood flows on the Grand, Conestogo, and Speed Rivers. 
The locations of the reservoirs are shown in Figure 2-1. 

Several aspects of the flow regulation contribute to non-homogenity in 
the sample of annual flood flows: 

a) available flood storage has increased as reservoirs were 
constructed. 

b) flood storage initially available for each event varies according 
to the weather conditions for months preceding the event. 

c) reservoir operating strategies have change over the years as 
experience and access to real-time data have increased. 

In order to carry out the frequency analysis of a homogeneous 
sample, it was thought important to estimate the natural flood flows 
(i.e. the flows which would have occurred had the reservoirs not been 
in place). The effects of the reservoirs could then be imposed in a 
consistent manner with the use of a reservoir operating model. The 
estimation of natural flows and the frequency analysis of natural 
flows will be discussed in this report. The effects of the reservoirs 
will be shown. The assumptions made and the modelling techniques used 
in determining the effects of the reservoirs will be discussed in 
Technical Report No. 38, "Reservoir Operations". 

2.3 Estimating Natural Flows 

The event which produces the annual maximum flow under natural 
conditions is not necessarily the same one that produces the observed 
maximum. Nor does a single event produce the maximum annual flow at 
each location of interest. The events that produced the annual 
maximum under natural conditions were selected by inspection of the 
natural daily flows (See Tech. Report No. 34). 

For each event selected, the following data were collected: 

a) hourly streamflow tables or stage charts for the Grand River at 
Gait and Brantford from Water Survey of Canada (WSC) or the Grand 
River Conservation Authority (GRCA); 
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b) hourly streamflow tables or stage charts for the Grand River 
below Shand Dam from WSC or GRCA, or operators' logs on gate and 
valve openings at the Shand Dam from GRCA; 

c) operators' logs on gate and valve openings at the Conestogo Dam 
from GRCA; 

d) lake level recording charts or operators' logs on lake levels at 
Shand and Conestogo Dams from GRCA; 

e) stage-storage and stage-discharge curves from Shand and Conestogo 
Dams from GRCA. 

The reservoir inflows were calculated based on the hydrologic 
continuity equation as follows: 

inflow = discharge + (change in lake storage) 

However, to produce a smooth hourly hydrograph from input data at 
irregular intervals, the inflows were calculated by: 

a) attributing the average inflow to the midpoint of the interval. 

b) connecting these midpoints when trend is continuing in one 
direction. 

c) connecting an extrapolation from each midpoint when trend changes 
during the interval . 

The volume of inflow is strictly preserved. 

The natural flow hydrograph for each event at Gait was estimated by 
replacing the portion attributed to reservoir discharge with the 
inflow. The steps involved were as follows: 

a) reservoir discharges from Shand and, if applicable, Conestogo Dams 
were routed, using the Muskingum routing technique, to Gait. 

b) the routed discharge hydrographs were subtracted from the observed 
hydrograph at Gait to produce a new hydrograph of local flows. 

c) the reservoir inflow hydrographs were routed to Gait and added to 
the local hydrograph. 

Muskingum K and X equal to 17 hours and 0.4 were used to route from 
both Shand Dam and Conestogo Dams to Gait. The values of these 
parameters were determined based on recent high flow events. 

The natural flow hydrogrpah for each event at Brantford was estimated 
in a similar manner by replacing the portion attributed to the 
observed flow from Gait by the natural flow at Gait. Muskingum K and 
X equal to 5 hours and 0.4 were used. 

Since 1974, natural hydrographs have been produced at all of these 
stations by the real-time flood forecasting model used by GRCA. The 
procedure is the same; however, streamflow data at many more gauges 
are included, thereby providing an even more reliable estimate. The 
natural flow hydrographs available by this means were used. 

Much of the required data for the 1950' s was not available. At that 
time, daily data was more generally used so operator's logs were not 
as complete. Some recording charts for major events were lost during 
follow-up studies. The hourly flow hydrographs could not be estimated 
using the above procedures. 



TABLE 2-2 - ANNUAL MAXIMUM INSTANTANEOUS NATURAL FLOWS 



Year 



1948 
1949 
1950 
1951 
1952 
1953 
1954 
1955 
1956 
1957 
1958 
1959 
1960 
1961 
1962 
1963 
1964 
1965 
1966 
1967 
1968 
1969 
1970 
1971 
1972 
1973 
1974 
1975 
1976 
1977 
1978 
1979 

1954 



f ows pr 



SHAND DAM 



~UdXe 



Mar 21 
Mar 28 
Apr 5 
Mar 31 
Apr 2 
May 2 
Mar 2 
Apr 3 
Apr 4 
Feb 27 
Mar 31 
Apr 8 
Apr 3 
Mar 28 
Mar 30 
Mar 27 
Apr 8 
Apr 12 
Mar 19 
Apr 3 
Mar 28 
Apr 18 
Apr 14 
Apr 12 
Apr 18 
Mar 12 
Mar 7 
Apr 19 
Mar 21 
Mar 14 
Apr 12 
Apr 14 



Flow 



25800 

7870 

12400 

10200 

8520 

4730 

9900 

11700 

13300 

4490 

3280 

8050 

11200 

3460 

10400 

9680 

4530 

5980 

3220 

8510 

5580 

5710 

8300 

8290 

14100 

7130 

14100 

15600 

14100 

13100 

8120 

12500 



C0NEST0G0 DAM 



Date 



Mar 20 
Mar 23 
Apr 4 
Mar 30 
Apr 1 
Mar 16 
Feb 16 
Apr 2 
Apr 4 
Feb 26 
Mar 29 



Apr 5 
Apr 3 
Feb 25 
Mar 30 
Mar 28 
Apr 8 
Feb 10 
Mar 18 
Apr 3 
Mar 21 
Apr 18 
Apr 14 
Apr 12 
Apr 17 
Mar 11 
May 17 
Apr 19 
Mar 21 
Mar 13 
Apr 11 
Apr 14 



Flow 



12600 
7900 
9410 
7110 
6750 
4630 
9890 
4960 
9650 
4620 

3500 

6350 

13000 

3650 

12800 

11700 

8340 

15600 

3130 

12000 

7440 

8280 

6720 

5250 

8590 

8990 

16650 

20800 

17200 

15320 

8460 

13600 



nted in italics are estimated 



GALT 



Uate 



Mar 20 
Mar 23 
Apr 5 
Mar 31 
Apr 2 
Mar 16 
Apr 17 
Mar 12 
Apr 5 
Feb 27 
Mar 29 
Apr 6 
Apr 4 
Feb 25 
Mar 30 
Mar 28 
Apr 8 
Apr 8 
Mar 19 
Apr 3 
Mar 21 
Apr 19 
Apr 15 
Apr 13 
Apr 19 
Mar 12 
May 17 
Apr 19 
Mar 22 
Mar 14 
Apr 12 
Apr 14 

Oct 16 



Flow 



57600 
30100 
35900 
30300 
21400 
13600 
23900 
17230 
31100 
12200 

7210 
32900 
46300 
10400 
29700 
28200 
15000 
26700 

9380 
37000 
14600 
23800 
17900 
21400 
32400 
28200 
55200 
58700 
42800 
40400 
26600 
36000 

59800 



based on dai ly 



BRANTFORD 



Hate" 



Mar 20 
Mar 23 
Apr 5 
Mar 31 
Apr 2 
May 3 
Mar 2 
Mar 12 
Apr 5 
Feb 27 
Mar 29 
Apr 6 
Apr 4 
Feb 25 
Mar 30 
Mar 28 
Apr 8 
Apr 8 
Mar 19 
Apr 4 
Mar 21 
Apr 19 
Apr 15 
Apr 12 
Apr 19 
Mar 12 
May 17 
Apr 20 
Mar 22 
Mar 14 
Apr 12 
Apr 15 

Oct 17 



mean 



Flow 



61100 
29800 
53100 
36200 
25500 
18600 
38300 
26000 
40500 
20300 

11300 
35700 
55500 
14100 
34600 
37000 
17000 
33600 
9320 
36500 
18400 
24300 
17600 
24300 
37900 
33800 
64700 
64000 
52200 
53200 
33100 
52700 

64700 



ows. 



GUELPH 



Date Flow 



Apr 5 

Mar 1 
Apr 1 
Mar 4 
Apr 17 
Mar 11 
Apr 5 
Apr 6 
Apr 7 
Apr 6 
Apr 4 
Feb 24 
Mar 30 
Mar 27 
Apr 8 
Apr 8 
Mar 19 
Apr 3 
Mar 20 



Apr 5 
Apr 9 
Apr 10 
Apr 15 
Mar 12 
May 17 
Apr 19 
Mar 21 
Mar 14 
Apr 12 
Apr 14 

Oct 16 



6720 
3780 

1870 
1660 
2870 
2420 
3670 
1930 

712 
3420 
3010 
1910 
2170 
3000 
1020 
3460 

920 
4810 
2300 
2750 
1640 
1700 
3730 
2930 
5850 
4480 
3500 
3700 
2280 
3340 

5740 



NEW HAMBURG 



Date Flow 



Mar 31 
Mar 22 

May 26 
Mar 2 
Apr 2 
Apr 4 
Feb 26 
Mar 28 
Apr 6 
Apr 4 
Feb 25 
Mar 30 
Mar 26 
Apr 7 
Feb 11 
Feb 12 



Apr 
Feb 
Apr 
Apr 



Apr 10 
Apr 13 
Mar 12 
Mar 5 
Apr 19 
Mar 21 
Mar 13 
Apr 11 
Apr 14 

Oct 16 



5570 

3630 

3630 

6270 

4530 

7650 

4060 

3170 

5990 

8820 

2470 

7010 

9220 

3750 

9540 

2230 

11800 

9880 

6950 

4440 

3680 

6290 

6800 

11000 

16000 

11400 

12800 

7030 

13700 

12000 



CANNING 



TJate" 



Mar 21 
Mar 23 
Apr 5 
Apr 1 
Mar 22 
Mar 17 
Feb 17 
Mar 12 
Apr 5 
Feb 27 
Mar 28 
Apr 3 
Apr 5 
Feb 23 
Mar 31 
Mar 37 
Apr 8 
Feb 11 
Feb 12 
Apr 4 
Feb 3 
Apr 6 
Apr 10 
Apr 4 
Apr 14 
Mar 13 
Mar 6 
Apr 20 
Mar 22 
Mar 14 
Apr 12 
Apr 15 

Oct 17 



Flow 



14900 

29800 

12600 

6720 

3970 

4590 

12900 

8940 

9420 

4850 

2360 

6510 

10200 

6500 

7320 

10700 

3610 

12500 

3540 

10400 

7040 

7010 

5100 

4280 

7350 

7160 

9710 

14800 

11100 

12800 

7900 

13500 

15100 



i 

en 

I 
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In order to obtain a continuous sample of annual maximum instantaneous 
flows, correlation equations were developed between maximum 
instantaneous and corresponding daily mean flows. The correlation 
equations and the plots of maximum instantaneous versus maximum daily 
flows are shown in Appendix A. 

The maximum instantaneous natural flows are summarized in Table 2-2. 

THe maximum instantaneous natural flows versus the corresponding 
observed flows are shown in Appendix B. 



3. FREQUENCY DISTRIBUTIONS AND ANALYSIS 

The frequency analysis was prepared using the Flood Damage Reduction 
Program Flood Frequency Analysis program (FDRPFFA) developed by 
Condie, Nix, and Boone, Inland Water Directorate, Environment Canada 
1977. 

The program gives the flood frequency regime as obtained using the 
Gumbel I, the Log-Normal, the Three Parameter Log-Normal, and the 
Log-Pearson Type III distributions using the maximum likelihood method 
of fitting. A moment fit is used when the maximum likelihood 
estimators do not have true solutions and is used for the Log-Pearson 
Type III distribution in addition to the maximum likelihood method. 
The asymptotic standard error of estimate is given as a percentage for 
selected return periods. 

The program was written in FORTRAN IV. For use on the IBM 5100 
desktop computer, it was converted to APL. 

A description of all distributions and a discussion of characteristics 
to be considered to determine reasonableness of fit are included in 
Appendix C. 

3.1 Flood Frequency Analysis 

Tables 3-1 to 3-14 show the flood frequency analysis using the 
four frequency distributions, as follows: 

Tables 

3-1 and 3-2 Natural Flow-Grand River at Shand 

3-3 and 3-4 Natural Flow-Conestogo River at Coenstogo Dam 

3-5 and 3-6 Natural Flow-Grand River at Gait 

3-7 and 3-8 Natural Flow-Grand River at Brantford 

3-9 and 3-10 Natural Flow-Speed River below Guelph 

3-11 and 3-12 Nith River at New Hamburg 

3-13 and 3-14 Nith River at Canning 

3.2 Selection of Flood Frequency Distribution 

In order to compare the fit of the Gumbel I, Log-Normal, and 
Three parameter Log-Normal distributions, the coefficients of skew and 
kurtosis of the sample were compared with the theoretical values of 
these coefficients. In the case of the Gumbel I distribution, the 
coefficients of skew and kurtosis of the sample should approach 1.14 
and 5.4 respectively. In the case of the Log-Normal distribution, the 



Table 3-1. Data Base and Plotting Positions 



NATURAL FLOW - GRAND RIVER AT SHAND DAM 
MAXIMUM INSTANTANEOUS FLQU (CUBIC FEET PER SECOND) 



194-3 - 1979 



YEAR 



DATA 



ORDERED 



RANK 



PROB 



RET. PERIOD 



19H-8 
19^+9 

1951 

1952 

T75T 

1954- 

1955 

T753" 

1957 

1958 

T93T 

1960 

1961 

1963 
1964- 



1966 
1967 
1963 
1969 
1970 



1971 
1972 
1973 
197i+ 
1975 
1976 
T97T 
1978 
1979 



25300 
7370 

124-00 

10200 

8520 

"T+73TT 

9900 

11700 



T33QO 
4-4-90 
3280 
3050 

11200 
3*160 

104-00 
9680 
1+530 
5990 
3220 
S510 



DuiJU 

5710 
3300 



8290 

14-100 

7130 



14-100 
15600 
14100 



lJTQTT 

8120 

12500 



25300 
15600 



14-100 
14-100 
14-100 
13JUU 
13100 
12500 



124-00 
11700 
11200 



104-00 

10200 

9900 

9480 
8520 
8510 
3300 
3290 
8120 
30SO 
7870 
7130 



5980 
5710 
5580 
4-730 
4-530 
4-4-90 



34-60 
3230 
3220 



o 

-> 

3 

""9" 
10 
11 

IT 
13 
14- 

TS" 
16 
17 

T5" 

19 
20 



23 
24- 



JO 

23 

19 



IT 
31 
32 



030 
061 



091 
121 
152 

212 



303 

333 



o6T 

394 

424 



4-85 
515 

3i+=r 

576 
606 

667 
697 



758 
733 

Tl& 

84-3 
879 

"WW 

939 
970 



33. 00 
16.500 



11 . 000 
3.250 
6.600 



b . i u u 
4-. 714 
4-. 125 



3,667 
3.300 
3. 000 



2.750 
2.533 
2,357 
2.200 
2.063 
1 .94-1 



1 .333 
1 .737 
1 .650 



1 . Zj < 1 

1 . 500 
1 . 4-35 



1 , 375 
1 , 320 

1 TiQ 



1 .179 
1 , 133 



1.100 
1 . 065 
1 , 031 




— 


NAIURAL FLOW GRAND RIVER Al SHAND MAM 
MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 
194B - 1979 

SAMPLE STATISTICS 




















— 


MIAN - 94 95 S.D. - 4 637. 9 




C.S. ■ 


1.3024 


C . K . *> 


6.0325 
3. 1229 


— 




SAMPLE STATISTICS (LOGS) 
MEAN - 9. 0446 S.D. - .4967 




C.S. - 


".2704 


C.K. - 


— 


SAMPLI MJN - 3220 SAMPLE MAX - 2580 Q 


N = 32 












PARAMETERS FOR LOUNORMAL M - 9.0446 


11 - 

S - 


/H50 
.4967 












PARAMETERS EUR IIIKEE PARAMETER LOGNORMAL 


A ■ 


"1824 


M = 9. 


2575 S ■ . 399t 


3. 1979 


-- - 


1 


STATISTICS OF LOG(X-A) 
MEAN - 9.25/5 S.D. ■ . JVVB 




C.S. = 


.0579 


C.K. - 


— 


PARAMETERS FOR LOO PEARSON III BY MOMENTS . A-= " 
PARAMETERS FOR LOG PEARSON IITBY MAXIMUM LIKFLIUO 
DISTK1DUI ION STATISTICS MEAN ■ 9.0446 S.D 


.0691 B«- 


51 .6098 


L0G(M)^12.6126 


M ■ 
LOG(M) -1 


300311 

2.2673 M - 


212632 


OS A ■ 

* .4091 


.0742 
C.S. 


B ■ 43.4179 
■ ".3035 


- 


















GOMDEL I LOGNORMAL 


1IIREE PARAMETER 




LOG PEARSON III 




LOGNORMAL 
RETURN FLOOD SI. ERROR FLOOD ST. ERROR FLOOD ST. ERROR 
PERIOD ESTIMATE PERCENT ESTIMATE PERCENT ESTIMATE PERCENT 

1 .005 1560 2090 2590 
1.050 3520 3/40 3610 


MAXIMUM 
FLOOD 
ESTIMATE 

2600 
364 


LIKELIHOOD 
ST. ERROR FLOOD 
PERCENT ESTIMATE 

2660 
36 1 


MOMENTS 
ST. ERROR 
PERCENT 








1 . 2511 5//0 5500 

2.00 G/4 04 70 

5.000 12/00 0.02 12900 10.20 

10.000 15400 9.37 16000 11.90 

20 . (I I /VO 9 . 84 1 920 13.50 

50.001) 21200 10.30 23500 15.50 


5660 
066 
1290 
15700 
104 
220 


9.3V 
10.90 
13.00 
16.00 


5660 
0600 
12900 
15600 
101 00 
21300 


5620 

06 70 

9.27 1290 

10.30 15700 

12.00 10400 

15.30 21(10 


15.60 
19.30 
24 .50 
32.70 








10(1.0011 2.5/00 10.60 26900 16.90 
200.000 26100 10.90 30500 10.20 
50 0.000 2V4 11.20 .554 1)0 19.90 


2TT700 
27500 
31300 


IB. 30 
20.6(1 
23.60 


23700 
26000 
29000 


19.20 24300 

21.20 2600(1 
25.50 30000 


39.70 
4 7.10 
57. 30 







J 
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Table 3-3. Data Base and Plotting Positions 



NATURAL FLOW - C0NE3T0G0 RIVER AT CONE3TOGQ DAM 
MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 
19 4 8 t9^r9 ■ 



YEAR 



1948 
19if9 

1950 



1951 
1952 
1955 



195** 

1955 
195 6 



1957 
1958 



1960 
1961 
19 6 2 



1963 
1964 

19 6 5 



1966 
1967 
19 6 3 - 



1969 
1970 
1971 



1972 
1973 
1974 



1975 
1976 
1977 



1973 
1979 



DATA 



12600 
7900 

— OJ 



44-e- 

7110 
6750 
'1 650 



9890 
4960 

9650 



4620 
3500 

&550 



13000 

3650 

12SQ0 



11700 
8340 

15 6 00 



3130 
12000 
-7'J i| Q 



8230 
6720 

525Q 



3590 
3990 



20300 
17200 
9- 



uji. 1 

3460 
13600 



ORDERED 



20800 
17200 

16550 



15600 
15320 



- J J O I 

13000 
12800 
12600 



12000 
11700 

9990 



9650 
9410 

8990 



8590 
8460 

95M3 



3280 
7900 
7'Ml 



7110 
6750 



6350 
5250 

4 960 



4630 
4620 
365 » 



3500 
3130 



RANK 



10 
11 
4-3- 



13 
14 



16 

17 

4^- 



19 
20 

-34r 



23 

-314- 



-.-I 



-32- 
2S 
29 

-5-9- 



31 
32 



PROD 



.030 
. 061 
t^4t 



.152 

i 132 



.212 



.303 
. 333 

644- 



I ■ >» ' 



.394 
.424 

.435 
.515 



' ■J 't o 

576 
606 



■ i OO 'O - 

. 667 
.697 



758 

733 

-34-3- 



.343 
.379 

.939 
.970 



RET. PERIOD 

33. 000 
16.500 



1 1 . 
3.250 
6.60Q 



4. 714 
4.125 



3.300 
3. 000 



2.533 
2.357 
2,200 



2.063 
1.941 
4 



1 .737 

1 . 650 



-i . •- ■ f r - 

1 .500 
1 .435 



X I ' J ' ' -> 

1 .320 
1 .269 

■J T*> **> 



1 . 179 
1 . 133 
44- 



1 . i U O ' 

1 . 065 

1 . 031 



-H- 



v Table 3-4. Flood Estimates 



NATURAL FLOW CONES TOGO RIVER A! CONESIUOO HAM 
MAXIMUM INST ANT AMEOU8 FLOW (CUBIC FEET PER SECOND) 
1V4U - iv/v 

SAMPLE STAT lb IICS 



MEAN - 9525 S.D. <* 4429.8 

SAMPLE STATISTICS (LOUS) 

MEAN - 9.0513 8.0'. " .4908 



C.S. - 



.6534 



C.K . - 



C.S. ■ 



.2831 



C.K. 



$.21 6B 



2. 744 6 



SAMPLE MIN - 3130 SAMPLE MAX - 20800 N ■ 32 

PARAMETERS FOR GUMDEL-I — -A — -.000285 U - 7H8U- 

PARAMETERS FOR LOONORMAL M - 9.0513 S - .4908 



PARAMETERS FOR (THREE PARAMETER LOONORMAL 

STATISTICS OF LOG(X-A) 

MtAN - 9.2973 8.0. « .3902 



1773 



M « 9.2573 



.3982 



C.S. 



. 1082 



C.K . ^ 



2.6381 



PARAMETERS FOR LOO PEARSON III DY MOMENTS A- .0695 D«* 49.9041 LOO(M )-12 . 5183 M m 27330 

PARAMETERS FOR LOG PEARSON III~D¥- MAXIMUM-LIKELIHOOD A = " . 1 387 —fr ■« 12.4843 LQG(M> -10.7826 

DISTRIBUTION STATISTICS MEAN - 9.0513 S.U - .4900 C.S. = ".5660 



48176 



COMPEL I 



LOONORMAL 



IHREE PARAMETER 
LOONORMAL 



III 



RF IORN 
PERIOD 

. 005 
. 050 

,2S0- 

.000 

5.000 

10.001) 

20 .000 

50 .000 

10U.000- 

20 0. 000 

50 0. 000 



FLOOD 
ESTIMATE 



16 30 

358(1 

— 5U1U 

8770 

120 

15400 

1 79 

21200 

-23400- 

26100 

2930 



ST. ERROR 
PERCENT 



FLOOD ST. ERROR FLOOD ST. ERROR FLOOD ST. ERROR FLOOD ST. ERROR 
ESTIMATE PERCENT ESTIMATE PERCENT ESTIMATE PERCENT ESTIMATE PERCENT 



LOO PEARSON 

MA XI M U M LIK E LI HOOD 

FLOOD ST. ERROR FLOOD 



MOMENTA 



2950 
380 

&64TL 



2650 
3670 



2500 
3550 
b/60 



2710 

366 



a. ii 

9.33 

9.80 
10.30 
10^60- 
10.90 
11 .20 



8530 
12V0 
16000 
1910 
23300 
3 6 700 
30 20 
35100 



10.10 
11 .70 
13.30 
15.30 
16. H 



18.00 

19. /0 



07 10 
12900 
15/00 
1O400 
21900 
4*4/00- 
2/500 
31200 



9.32 
10.90 
12.90 
15.90 

-TO .20 



20.50 
23.4 



8930 
13000 
15400 
17500 
20000 
21700 
2330 
25300 



8.51 

8.86 

9.88 

12.00 

-E4-HF©- 



16.20 
19.30 



5690- 
8730 
13000 
15700 
18300 
21600 

2650 
29700 



15.60 
19.20 
24 . 4 
32.60 
-39 .6 



47.0 
57.20 
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Table 3-5. Data Base and Plotting Positions 



NATURAL FLOW - GRAND RIVER AT GALT 
MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 

1948 - 1979 



YEAR 



DATA 



ORDERED 



RANK 



1948 
1949 



57600 
30100 
35900 
30300 
21400 



PROB. RET. PERIOD 



53700 
57600 



. 030 
, 061 



J<3 
16 



00 
500 



195 
1951 
1952 

T7ST 
1954 
1955 
1956 
1957 
1958 

1757" 
1960 
1961 
1962 
1963 
1964 
lVob 
1966 
1967 
1963 
1969 
1970 
1971 
1972 
1973 

"IV ,'i+ 



55200 
46300 
42800 
40400 
37000 
36000 
3590(1 
32900 
32400 



, 091 
.121 
. 152 
"7TZ7 
. 212 
51*2 



11,000 
8.250 
6. 600 



1975 
1976 
TTTT 
1973 
1979 



lodUU 
23900 
17230 



31100 
12200 
7210 
3290(1 
46300 
10400 



29700 
28200 
15000 
26/0 
9330 
37000 
14600 
23800 
17900 
21400 
32400 
28200 



58700 
42300 
4-U4UG 
26600 
36000 



31100 
30300 
30100 
2970(5 
23200 
2S200 
24700 
26600 
23900 



23300 
21400 
21400 



17900 
17230 
15000 



14600" 
13600 
12200 
10400 
9380 
7210 



9 

10 
11 



12 
13 
14 



15 
16 
17 



13" 
19 
20 



22 
23 



24 
25 
26 



29 



31 
32 



.273 
.303 

.333 



.364 
.394 
.424 



.455 
.485 



jij 



34T 
.576 
. 606 



. 606 
. 667 
.697 



.753 
■ 783 
.313 
. 843 
,879 



. 939 
. 970 



5.500 
4. 714 
4. 125 



3.667 
3.300 
3. 000 



"7*.\ 



u0 
2.533 
2.357 



2. 200 
2.063 
1 .941 
1 .333 

1.737 
1 . 650 



371 
500 
435 



320 
269 



1 . 222 
1 . 179 

1 . 133 



1.10 
1 . 065 
1 . 031 



> 



4-i- 



Table 3-6. Mood Estimates 



NATURAL FLIJU GRAND RIVER A I BALI 

MAXIMUM INSTANTANEOUS FLOW (CUBIC FEE] PER SECOND) 

" 1948 1979 

SAMPLE STATISTICS 



< 



MEAN - 28566 S.B. - 

SAMPLE STATISTICS (LOGS) 
MEAN - 10.1365 S.U. - 



13638. / 



C.S. - 



. 5836 



C.K. - 



>301 



C.S. = 



.5431 



C.K. - 



3.2720 



J. 14 76 



SAMPLE MiN - /210 SAMPLE MAX - 58700 N - 32 

PARAMETERS FOR GUMBEL I A"- .000090 U - 52~224~ 

PARAMETERS FOR LOGNORMAL M - 10.1365 S - .5301 



PARAMETERS FOR THREE PARAMETER LOGNORMAL 

STATISTICS OF LOG(X--A)' 
MtAN - 10.6032 S.U. ■ .30 0'* 



A - 



16980 



M ■ 10.6832 



S ■ .3004 



C.S. * 



.0681 



C.K. ■ 



2.7912 



PARAMETERS FOR LOG PEARSON III BY MOMENTS A- .1439 6= 13.5622 LOG(M ) =12 . 0886 

PARAMETERS FOR LOG PEARSON III BY MAXIMUM LIKELIHOOD A = ".25'OT 6 • 

DISTRIBUTION STATISTICS MEAN - 10.1365 S.D ■ .5334 C.S. - .8251 



M ■ 177825 

5.8757 LOG(M) =11. 4 295 



91998 



GUMBEL I 



LOGNORMAL 



THREE PARAMETER 



LOG PEARSON III 



Rb TURN 
PERIOD 

1 .0 05 
1_, 050 

i :25o 

2.000 

5 . 

10.000 

20 . 00 

50 . 

100 . 00(7 

20 . 00 
50 .0011 



FLOOD 

EG I IMA IE 



35 V 
9790 



16900 
26 50 
.590 
'I /'I 
55400 
6580 

TSEM 

O i n o o 

9 I 6 



ST .ERROR 
PERCENT 



9.13 

9 . 6'i 
III 10 
10.60 



T(T8T 

11.10 
I 1 '10 



FLOOD ! 
ESTIMATE 

8010 
10608 



LOGNORMAL MAXIMUM LIKELIHOOD MOMENTS 

I. ERROR FLOOD ST. ERROR FLOOD ST. ERROR FLOOD ST. ERROR 
PERCENT ESTIMATE PERCENT ESTIMATE PERCENT ESTIMATE PERCENT 



1620 
25200 
.5950 
•l 98 
60400 

74900 

86600 

989 

I 16 II 



10.90 
12.60 
14.1 
16.50 
18.80 
19.50 
21 .30 



5/80 
96.50 



16900 
2660 
39200 
47100 
54500 
_6380 0_ 
70 .'0 
//600 
86600 



8.86 
10.00 
11 .70 

14 . :»o 



16.20 
18.20 
20.80 



6150 
94 50 



6750 

98 1 



16/00 
27100 
3980 
4690 
5280 
5920 



8 . .U, 
8.18 
8.51 

9 . 64 



1650 
2650 
39/0 
4/900 
5520 
64 000 



17,20 
22.0 
28.90 
4 0.0 1) 



63300 
6690 
7 1 



10.90 
12.40 
14 .50 



70 100 
7580 
829 



49.4 
59.30 
/.5. 10 
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Table 3-7. Data Base and Plotting Positions 



NATURAL FLOW - GRAND RIVER AT BRANTFORD 
MAXIMUM INSTANTANEOUS FLQU (CUBIC FEET PER SECOND 
1948 - 1979 



YEAR DATA ORDERED RANK PROS. RET. PERIOD 



l?i+8 61100 61+700 1 .030 33.000 

1949 29300 64000 2 .061 16.500 

•±2*& 52UU2 611QQ 3 (L9J n nnn 

1951 36200 55500 4 ,121 3 ->50 

1952 2550 S320 5 .152 6.' 60 
-i253 13AM 531QQ 6 %B2 * ^nn 

1954- 38300 52700 7 .212 

1955 26000 52200 8 .242 

-1254 MUM Mnunn 9 ' yr 



4,714 
4.125 

c.6 7. 



1957 20300 38300 10 .303 3.300 

J2I i 13 °° 3790 ° U ' 333 3.000 
■1212 257JIQ 37Q0fl 12 Z£k 2 7«rn 

1960 55500 36500 13 .394 2 538 

1961 14100 36200 14 ,4^4 2*357 
±242 ZMM 337 fin 15 ^55 ?'?nn 

1963 37000 34600 16 .485 2 063 

1964 17000 33800 17 .515 l.*94i 
4?4? 23AM 23AM 13 5k5 \ ' 335 



1966 9320 33100 19 .576 1 73', 

\lil ?i 5 ? n ° 2980 ° 20 - 606 i»«fl 

-^SAS IffltiLO 2ASLM 3J 4^ 1|5?1 



19o9 24300 25500 22 .667 TToo 

197 ° 17600 24300 23 6?7 1*433 

-1221 ?4^nn ^nn nu. %£ 



1972 37900 20300 25 .758 ~ 1.320 
197o 008OO 18600 26 .783 l ^aq 

■±22± t&JM lflunn n 7 aiQ « ' ;g; 

!ot? S!!; 1760 ° 2S -8*8 1.179 

1976 52200 17000 29 . 87* 1 iiq 

-1222 522JLQ ULLM ^n qhq j'jnn 

1973 33100 U300 31 .939 1 065 
1979 S2700 9320 32 .970 1 031 




f NATURAL FLQM GRAND RIVER Al BRAN 1 FORD 

MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 
1VUQ 1272 

SAMPLE STATISTICS 








- 










MEAN ■ 346V4 8.D. ■ 13685.3 
SAMPLE STATISTICS (LUGS) 




c.s. ■ 


. 3565 


C.K. - 


2.5543 
3.0631 




MIAN 10.341b Ji . 1 1 . * .3866 




c.s. - 


'.5525 


C.K . ■ 


SAMPLE MIN - 2321 SAMPLE MAX - 64 70 
PARAMETERS FOR UUMfctL I A - .000076 II - 


N - 32 
27261 






- 




PARAMMERS FOR LOGNORMAL M - 10.3m5 


S - 


.5066 






PARAMETERS FOR THREE PARAMETER LOtiNORMAL 

siai isrics OF i.oo<x-a> 


A = 


"27183 


M - 11. 


0015 S = .2553 


} 
2.5120 


i 
£ 


MtAN - 11.0015 S.D. « .2559 




C.S. ■ 


". 0764 


C.K. ■ 


PARAMEIFRS FOR LOG PEARSON III BY MOMENTS 
PARAMETERS LOR LUli. PtARSUN 111 BY MAXIMUM 
DISTRIBUTION STATISTICS MEAN - 10.3415 


A- .1400 B- 
.IKt IIMOOB A m 


13. 1029 
" . 3647 


L0G<M>=12. 1753 
B b 2 3117 


M - 193933 
LOG(M) ^11 .1845 M - 


/20 V 


S.B = .5545 


C.S. 


«"t .3154 












MOMENTS 




CUMULI. I LOGNORMAL 


IMREE PARAMETER 

i nnunoMAt 


MAXIMUM 


LOG PEARSON III 
i tkfi IMOOB 


RETURN FLOOD ST. ERROR FLOOD ST. ERROR 
PERIOD ESTIMATE PERCENT ESTIMATE PERCENT 


FLOOD 
ESTIMATE 

7260 
1220 
21200 


ST. ERROR 
PERCENT 


FLOOB 
ESTIMATE 


ST. ERROR FLOOB 
PERCENT ESTIMATE 


ST. ERROR 
PERCENT 




1.0 05 530 10 300 
1.050 12600 13500 
1.250 21000 20 20 




6220 
10 700 
209110 


8760 
1250 
20600 


2.00 32100 31000 

5.000 4/000 0.92 4/500 10.40 

lo.ooo 56Voo y.46 59300 12.10 

20.00 664 9.92 /J 30 13,70 

50.000 /U600 10.40 G/600 15. GO 

100. Olid B7800 10. /0 101000 17-20 


3 20 

4 7280 

56100 
64 20 
/4 20 
G1600 


0.31 

_ 9.30 

10. B0 

13.20 

is. an 


31800 
4V20 
550 
60500 
64 70 
6/000 


3250 
6.69 4 70 
5.85 57100 
5.15 654 
4.41 75100 
4 . ' OVMMM) 


16.80 
21.60 
20.40 
39.38 

40 60 




20 0.000 97000 1100 114 000 1G.60 
500.000 109000 11.20 133000 20.30 


GOG 
90108 


16.90 
IV. 30 


6G6 
70 000 


3.80 GG30 
3. /6 V6100 


58 . 4 

72. 





-16- 



Table 


3-9. Data Base 


and Plotting Positions 








NATURAL FLOW - SPEED 
MAXIMUM INSTANTANEOUS FLOW 
- - 1950__tl 

DATA ORDERED 


RIVER BELOU GUELPH 
(CUBIC FEET PER SECOND) 
1979 


YEAR 


RANK 


PROS. 


RET. PERIOD 


1950 
1951 
1952 


6720 
3780 

... 1370 

1660 
2870 
2420 


6720 
5850 

«+aiji 


1 
2 

3 


. 032 
.065 
.097 


31 . 000 
15.500 
10 . 333 


1953 
195«+ 
1?55 


i+i+80 
3780 
3770 


4 
5 


.129 
.161 
.194 


7.750 
6.200 
5. 167 


1956 
1957 
1958 


3770 
1*60 

712 __ 

3m o 

3010 
li+.-?n 


3730 
3700 
3500 


7 
8 
9 


.226 
.258 
.290 


4.429 
3.875 
3.444 


1959 
1960 

1961 


3U60 

3*10 

331+0 


10 
11 
12 


.323 
.355 
.387 


3. 100 
2.813 
2.583 


1962 
1963 
1?64„ 


2170 
3000 
1020 


3010 
3000 
293Q 


13 
14 
15 


.419 
.452 
.481+ 


2.385 
2.214 
2. 067 


1965 
1966 

1?67 


3460 
920 

i+SiQ 

2300 
2750 
1640 
1700 
3730 
2?30 


2870 
2750 
24-20 


16 
17 
18 


.516 
.548 

.581 


1.938 
1.824 
1 . 7"^ n 


1968 
1969 
1970 


2300 
2280 
- .2170 
1870 
1700 
1660 


19 
20 
21 

23 
24 


.613 
.645 
.677 


1.632 
1 .55 
1 . 476 


1971 
1972 
1973 


.710 
.7»+2 
,774 


1 .409 
1 .343 
1 .292 


1971+ 
1975 
1976 
1977 
1978 
1979 


5850 
4-1+80 
_ 3500 
3700 
2230 


1640 
1460 

. 1430 

1020 
920 
71? 


25 
26 

27 

28 
29 
30 


. 306 
.839 
.871 
,903 
.935 
,968 


1 .240 
1 . 192 

1 , 143 

1 . 107 " 
1 . 069 
1 . 033 





Idble 3-10. Mood Estimates 




-> 




NA1URAL FLOW - SPEED RIVER BELOW GIJELPM 
MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 

itai - j 9/9 ^ ... 

SAMPLE STATISTICS 








MEAN * 2B90 S.D. - 1403.0 C.S. - 

SAMPLE STATISTICS (LOOS) 
MIAN - /.U444 S.D. - .5326 C.S. ■ 


7812 C.K. - 4.1762 
.5602 C.K. = 3.4246 






SAMPLE MIN = 712 SAMPLE MAX ■ 6720 N - 30 
PARAME TERS FOR GUMDEL I A - .000086 U • 2247 
PARAMETERS FOR LOGNORMAL M - 7.0444 S ■ .5326 


.... 






PARAMETERS FOR THREE PARAMETER LOGNORMAL A ■ "1534 

STATISTICS OF LOG(X-A) 

MEAN ••= U.3471 S.D. = .3151 C.S. ■ 


M = 8.3471 S - .3151 
".0505 C.K. = 3.0753 


• 




PARAMETERS FOR LOG PEARSON III BY MOMENTS A= ".1492 B= 
PARAMETERS FOR LOG PEARSON IJI BY MAXIMUM MKELIHQPP CL - 


12.7466 LOG<M)= 9.7460 M * 17085 

".W24 p = 9.3584 LOG(M) - 9.4575 M ■ 


1280 3 


DISTRIBUTION STATISTICS MEAN - 7.8444 S.D * .5273 


C.S. - ".6538 




. 




------- 


DUMBEL I LOGNORMAL THREE PARAMETER 

LOGNORMAL 


LOG PEARSON III 
MAXIMUM LIKELIHOOD MOMENTS 




RETURN FLOOD ST. ERROR FLOOD ST. ERROR FLOOD ST. ERROR 
PERIOD ESTIMATE PERCENT ESTIMATE PERCEN1 ESTIMA1E PERCENT 


FLOOD ST. ERROR FLOOD ST. ERROR 
ESTIMATE PERCENT ESTIMATE PERCENT 

664 674 

984 985 

1680 1660 




1.0 05 364 805 598 
1.050 990 1060 978 
1.250 1711 1630 WOO 




2.000 2660 2550 2680 

5.000 3941 9.43 4 000 11.30 39/0 9.29 

10.000 4790 9.95 5050 13.10 4/80 10.50 

20.000 5600 10.40 6130 14.90 5551 '12.30 

50.000 6650 10.90 7610 17.10 6520 15.00 

10 0.000 74N1 11.20 8810 18.70 7240 17.10 


2700 2680 
4010 9.15 4030 17.80 
4 780 9.34 4850 22.90 
5450 10.20 5590 30.20 
6230 12.10 6460 41.80 
6750 14.10 70 70 51.70 






200.000 8230 11.50 10100 20.20 7960 19.20 
500.000 9260 11.70 11800 22.10 B921 22.00 


7230 16.20 764 62.20 
7810 19.20 834 76.80 


,j 
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Table 


3-11. 


Data Base and 


Plotting Positions 










NITH RIVER AT 
MAXIMUM INSTANTANEOUS FLOW 

1 o«ri _ 


NEW HAMBURG 
(CUBIC FEET 

1°7° 


PER 


SECOND) 


YEAR 




DATA 




ORDERED 


RANK 




PROB 


. RET. PERIOD 


1951 
1952 
I«51 




5570 
3630 




16000 

13700 
1 ?qnn 


1 

m 
*• 

3 




. 033 
.067 

. 1 nn 


30. 000 
15. 000 
1 Q , n n 


1954 
1955 
1 oe ?4 




6270 
4530 




11800 

114-00 

1 1 nnn 


4 

5 




.133 

.167 

onn 


7.500 
6. 000 

? , n n o 


1957 
1958 
1959 




4060 
3170 




9880 

9540 
goon 


7 
8 




233 

.267 

7nn 


4.286 
3.750 
3 .77^ 


1960 
1961 
196° 




8820 
2470 

791.9 




8820 

7650 
-7mn 


10 
11 
1? 




.333 
367 
an n 


3. 000 
2 7 n 7 


1963 
1964 
19*5 




9220 

3750 

9540 




7010 
6950 
&SO&. 


13 
14 

15- 




433 
.467 

..5&Q. 


2.308 

2.143 
-> nnn 



1966 
1967 
1948 



1969 
1970 
1771 



1972 
1973 
197H 



1975 
1976 
1977 



2230 
11800 

9380 



6950 
4440 
3630 



6290 
6300 
,10 



16000 
11400 

12300, 



6290 
6270 



5570 
4530 
■a li 4 



16 
17 
-12- 
19 
20 



567 



600 

667 

7nn 



4060 
3750 
34843 



2 3 



767 
■&Q-0- 



3630 
3630 

3170 



26 



333 
867 

'?nn 



1.375 
1.765 

1 , 6 A 7 



1 .579 
1.500 

1 ,479 



1 . 364 
1 ,304 
i 2SJ3 



1 .200 
1 . 134 

1 111 



1978 
1979 



7030 
13700 



2470 
2230 



28 
29 



933 
967 



1 . 071 
1 , 034 



Table 3-12. Hood Estimates 



- 


IHIII KlVlk Al III U MAMhUKli 

MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SET 
1231 = 13.73 - 

SAMPLE STATISTICS 


:ond) 






\ 








— 


MEAN * 7218 S.H. - 3623.6 
HAM I'll SIANSIICH (LOGS) 


c.s. - 


.6730 


C.K. - 3.0 701 




MEAN - 11.7567 S . U . ■ .526V 


c.s. - 


". 1903 


C.K. = 2.5461 


— 


SAMPLE MEN 2230 SAMPLE MAX - 16000 N * 29 
PAKAMb It-'MK FOR GIIMHEI. I A - . HIHHr./ II - KSSO 








PARAMETERS FOR LOGNORMAL M ■ 0.7567 


S = .5269 






— 


PARAMETERS FOR THREE PARAMEIER LUGNORMAL 

siai ihi res in i nii< * a) 


A - "170 


M - 8. 


7068 S =■ .51 12 




MIAN - IJ.7M68 9.0. ■ .5112 


C.S. - 


-, i63y 


C.K. = 2.5246 


— 


PARAMETERS FOR LOS PEARSON III BY MOMENTS 
PARAMETERS EOH LOO PEARSON 111 BY MAXIMUM L 
DISTRITJUTION STATISTICS MEAN - 8.7567 


A^ .0501 B^ 

ik! i.iiinnii a 1 


110.4524 
.1413 


L00(M)-14.2942 M - 1613982 

B ■ 13.0621) I.OO(M) -18.7ISS M 


1 
45050 f 


S.D " .5261 


C.S. 


- ".5372 














GUMDEl. 1 l.OONORMAL 


IIIREE PARAMETER 
1 llliWUHMAI. 


MAXIMUM 


LOG PEARSON 111 

I.IKFI lllflllli MIIMFNTS 


- 


RL IIJRN FLOOD ST. ERROR H.UUH ST. ERROR FLOOD ST. ERROR 
PERIOD ESIIMAIE PERCENT ESTIMATE PERCENT ESTIMATE PERCENT 


FLOOD 
ESTIMATE 


ST. ERROR FLOOD ST. ERROR 
PERCENT ESTIMATE PERCENT 




1.0 05 OU4 2030 
1 .050 24 4 2670 
1 :>!-.« ii'MII hiuia 


IVVO 
2650 
41) VII 


1720 

24V0 
uun 


1910 
260 
4 1 




2.00 65V 6 350 

5.00(1 V/60 V.61 WOO 11.40 

10.000 11V00 10.10 12500 13.20 

20.0 0(1 13900 10 .60 15100 15.00 

50.000 16580 11.10 18/00 17.20 

i on non ni4 no nun .11 /.no i«,i'MI 


6300 

WOO 11.20 
12400 13.40 
15000 .16.00 
1050 19.60 
.'1300 22.511 


6660 
W70 
12000 
13000 
16000 
1 /Mill 


646 
9.6V VV4 16.60 
10.20 1230 20.10 
11.40 14700 25.20 
14.00 17700 33. HO 
16 Ul 2111 Oil 4 11 40 






200.000 20400 11.70 24/00 20.30 
500.000 2.IOO0 11. 90 29890 22.20 


24300 25.20 
2114 20.00 


19000 
20O0 


IV. 00 22500 4/. 70 
22.60 2570 57.90 
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Table 3-13. Data Base and Plotting Positions 



N1IH KIVLN Ai CANNING 



MAXIMUM INSTANTANEOUS FLOW (CUBIC FEET PER SECOND) 

1943 - 197? 



YEAR 



lV4ti 
191+9 
1930 
1951 
1952 
1953 
TTS*- 
1955 
1956 

195 r 

1958 
1959 

I 9 6 0' 



1961 
1962 
1965 - 
1964 
1965 
1V66 
1967 
1968 
" 1 9 6? 
1970 
1971 
TTT2 
1973 
1974 
1975 
1976 
1977 
T7T3" 



19" 



DATA 



1M-VUU 
9910 
12600 
— 6720 
3970 
1+590 
1290Q 
891+0 
9420 
1 4350 
2360 
6510 
102 
6500 
7320 
1030 
3610 
12500 
334U 
10400 
701+0 
— " 01 



5100 

4280 

7 3JU 

7160 

9710 

1430 -0- 

11100 

12800 

,-VUU 

13500 



ORDERED 



m-yuu 

11+800 
13500 
12 9 00 
12800 
12600 
12O0Q - 
11100 
1 OM-0 
10300 



10200 
9910 

9 710 - 



9420 
891+0 
7 900 
7350 
7320 
7 160 
701+0 
7010 
5 7 2 
6510 
6500 
31 JO 
4350 
4590 
i+230 
3970 
3610 
SZTtT 
2360 



RANK 



9 

-rt- 
11 
12 

-1-3- 

14 

15 

-w 

17 

18 

-rr- 



21 

■■ ■* m> 



24 



26 

27 

-»■ 

29 



"3r 

32 



-mnr 

061 
091 

-tei- 



152 
132 

q 4 r>. 



273 
-2Ffr3- 



.333 

.361+ 
, 39 1 4 



.455 

7*83- 

.515 

.545 

r3T3- 

.606 

636 

,667 

697 

727 

-7^3- 

738 

818 

843 



.37? 
.90? 
■ -7 5 y 
.970 



PROB. RET, PERIOD 



33 . U (J U 

16.30 

11. 000 

-■ 8.25 

6.600 

5.500 

4.714 



4. 125 
3 . 667 
3.3QQ 
3. 000 
2.750 

-1 rj-SQ. 



2.357 
2,200 
2. 063 
1.941 
1 .333 
1 .737 
1 .650 
1.571 
1 .GOO- 



1 .433 

1.373 
1 . 32 
1 .26? 
1 . '" , '" ,0 

* 1 "iq 



1 . 133 
1 . 100 

I . J 63 
1 , 031 



> 



Table J-14. Flood Estimates 



nttm rtvtr-at canning 



MAX 1 HUM lNiilANIANEOIJB FLOW (CUBIC FEET PER SECOND) 
1948 - 197V 



SAMPLE STATISTICS 
MFAN~= 8«ft« 



S. D . = 



3511 .3 



c .a. - 



.-taa^ 



C .K. 



2 .3? 25 



SAMPLE STATISTICS (LOGS) 

MIAN ■ 8.9424 S. Ii. - 

^5AMPLF~MTN = " 233U— 



.4704 



c.s. = 
tt-« — 32— 



.6045 



C.K. 



3.0630 



PARAMETERS FOR GOMBEL I 



.000328 



U - 



6740 



S = 

-A—=^ 



.4704 
" 143 4 



H = 1 . 2 



-S~= ■ 1552 



PARAMETERS FOR LOGNORMAL M - 0.9424 

^ARAMFTFRS-FOR T HR E E P ARAMETE R L OGNORtlAtT- 

S1AIISIICS OF Lf)G<X A) 
MtAN ■ 10.020 S.D. ■ .1552 C.S. ■ ".0529 C.K. ■ 2.3659 

PARAMrTERS-rOR-t-OO-ptARaON—ltt— tn—flOMtrNfS A- ~. 1422 B- tO . 9 47 L OO ( M ) - 1 ■ H » B » M~- 36276- 

PARAMETERS FOR LOG PEARSON III bY MAXIMUM LIKELIHOOD A - .4195 B ■ 1.1744 LOG(M> = 9.4351 

DISTRIBUTION STATISTICS MEAN- B.9424 S.D- .4546 C.S. -"1.8455 



M - 



12520 



I.IIMI'l I. 1 



LOGNORMAL 



— RETURN 
PERIOD 

1 . 005 
1 .050 
1 . 250 

2T000 

5.000 

10.000 

20 . II 

50.000 

10 0.00 

zo ot onn 

500.000 



FLOOD" 
ESTIMATE 



1660 

3350 

52 V 

-7860 

1 130 

13600 

15800 

10600 

20 /0 

22?00~ 

2570 



8.58 

9. 15 

9.64 

10.20 

10.50 

TO. 80 

11.10 



2760 
3530 
5150 



7630 

114 
14 
16600 
20100 
22000 
2570 0- 
29/0 



9 
11 
12 

14 

16 



60 
20 
8 
70 




1 7 . 30 
10. 90 



THREE PARAMETER 
LOGNORMAL 

:OOTi S T . ERROR 

ESTIMATE PERCENT 



LOG PEARSON III 
MAXIMUM LIKELIHOOD 

-FtOOD ST . E RR OR — P tOO B 

ESTIMATE PERCENT ESTIMATE 



MOMENTS 
3 T . ERROR 
PERCENT 



—1750 

3110 
54 20 
— 817U 

11300 

1310 

IN /0 

1660 • 

1 790 
T7^00 

20800 



7.25 
7.92 
9.12 
11.00 
12.60 
14. 1 
16. 20 



2330 
3280 
5260 



OH 
00 
00 
00 


-otr 
oo 



— 8*K*e- 

114 

13400 

15200 

1720 

18500 

TYtttttt- 

21300 



16.40 
21 .20 
28 . 1 
39. 10 
48.50 
-58-r5fr- 
72.40 
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coefficients of skew and kurtosis of the logarithmic transform with 
the addition lower boundary parameter should also approach 0.0 and 
3.0. 

The coefficients of skew and kurtosis for each station compare to the 
theoretical values as shown in Table 3.15. 

It can be seen form the comparison of the coefficients of skew and 
kurtosis that of the Gumbel I, Log-Normal and Three-Parameter 
Log-Normal distributions, the Three-Parameter Log-Normal provides the 
best fit. 

In all six cases, the fitting of the Log-Pearson Type III distribution 
produced an upper bound. The upper bound for both the moment and 
maximum likelihood fit at each station are shown below: 

Table 3-16 - Upper Bound of Log Pearson Type III Distribution 



Station 



Upper Bound by Ratio to Upper Bound Ratio To 
Max. Likelihood Sample Max. Est. By moments Sample Max. 



Shand Dam 


212632 


8.2 


300311 


12. 


Conestogo Dam 


48176 


2.3 


273300 


13. 


Gait 


91998 


1.6 


177825 


3.0 


Brantford 


72009 


1.1 


193933 


3.0 



Guelph 

New Hamburg 

Canning 



12803 



45050 



12520 



1.9 



2.8 



0.8 



17085 



1613982 



36276 



2.5 



100. 



2.4 



In the case of the maximum likelihood method of fit, the upper bound 
at five of the six stations sets an unreasonable restriction on the 
ordinate and produces much lower flows for low frequency events than 
do the other distributions. For Canning, the upper bound is lower 
than the sample maximum and the distribution cannot be used. 

In the case of the moment method of fit, the upper bounds are 
general ly much higher and have only marginal physical significance. 
The resulting flood estimates are similar to those produced by the 
Three Parameter Log-Normal distribution. It is felt that the 
Log-Pearson Type III with parameters estimated by moments provides a 
reasonable fit. 



A comparison of the fit of the Three Parameter Log-Normal and the 
Log-Pearson Type III distribution moment solution is shown for each 
station in Appendix D. 

A conclusion of a recent study done by Maclaren Engineering Planners & 
Scientists Inc., for The Ontario Ministry of Natural Resources, 
October, 1981 (5) states: 



TABLE 3-15 Comparison of Coefficients of Skew and Kurtosis 



Station Coefficient of: GUMBEL I 

Sample Theoretical 



LOG-NORMAL 3-PARAMETER LOG-NORMAL 
Sample Theoretical Sample Theoretical 



Shand Dam 


skew 


1.3024 


1.14 


-.2784 


0.0 


-.0579 


0.0 




kurtosis 


6.8325 


5.4 


3.1229 


3.0 


3.1979 


3.0 


Conestogo Dam 


skew 


.6534 


1.14 


-.2831 


0.0 


-.1082 


0.0 




kurtosis 


3.2168 


5.4 


2.7446 


3.0 


2.6381 


3.0 


Gait 


skew 


.5836 


1.14 


-.5431 


0.0 


-.0681 


0.0 




kurtosis 


3.2720 


5.4 


3.1476 


3.0 


2.7912 


3.0 


Brantford 


skew 


.3565 


1.14 


-.5525 


0.0 


-.0764 


0.0 




kurtosis 


2.5543 


5.4 


3.0631 


3.0 


2.5128 


3.0 


Guelph 


skew 


.7812 


1.14 


-.5602 


0.0 


-.0505 


0.0 




kurtosis 


4.1762 


5.4 


3.4246 


3.0 


3.0753 


3.0 


New Hamburg 


skew 


.6730 


1.14 


-.1903 


0.0 


-.1639 


0.0 




kurtosis 


3.0781 


5.4 


2.5461 


3.0 


2.5246 


3.0 


Cann ing 


skew 


.1865 


1.14 


-.6045 


0.0 


-.0529 


0.0 




kurtosis 


2.3725 


5.4 


3.0630 


3.0 


2.3659 


3.0 



t 

CO 

I 
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"Estimates of peak flow using the Log-Pearson type III Distribution 
are somewhat more sensitive to variations in log skew than are 
estimates from the Three Parameter Log-Normal Distributions for 
similar variations in real skew. Generally speaking, the Three 
Parameter Log-Normal Distribution was found to be marginally superior 
compared to the Log-Pearson Type III Distribution". 

This report makes the following recommendation: 

"The Three Parameter Log-Normal Distribution is suggested for flood 
frequency analyses in the Province. If the coefficient of skew of the 
untransformed series is negative or the Three parameter Log-Normal 
distribution cannot otherwise be fitted, then the Log-Pearson Type III 
distribution should be used". 

In view of the comparison of the distributions and the recommendation 
regarding the distribution to be used for future studies in Ontairo, 
the Three-Parameter Log Normal distribution was selected for use in 
this study. 

3.3 Period of Analysis 

In general, the period of record from 1948-1979 was used for the flood 
frequency analysis. (The period of record for Guelph and New Hamburg 
started with 1950 and 1951 respectively as the gauges were installed 
in those years). A common period of analysis at all stations was 
selected principally to aid in the development of flood frequency 
relationships for flood prone areas along the river 
(See 1.0 Background). 

For the Grand River at Gait, only the last 32 years of 66 years of 
record were included in the analysis. A comparison of frequency 
curves developed from 32 years of record (1948-1979) and from 66 years 
of record (1914-1979) is shown in Appendix E. The shorter term 
analysis resulted in 12-14% higher flow estimates than the longer term 
analysis for the same return period. The frequency curve for the 
shorter term analysis lies within the upper 95% confidence interval of 
the frequency curve for the longer term analysis. 

4.0 RESULTS 

Flood frequency analysis was performed on maximum instantaneous spring 
freshet flows at the following long term gauge stations in the Grand 
River Basin: 

Grand River at Shand Dam 1948-1979 

Conestogo River at Conestogo Dam 1948-1979 

Grand River at Gait 1948-1979 

Grand River at Brantford 1948-1979 

Speed River below Guelph 1950-1979 

Nith River at New Hamburg 1951-1979 

Nith River at Canning 1948-1979 

The effects of regulation by the Shand, Conestogo, and Guelph Dams on 
the Grand, Conestogo, and Speed River stations were removed. 

Where maximum instantaneous flows were not available, they were 
estimated based on the daily mean flow. 
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Analysis was performed using the Gumbel I, Log-Normal, Three Parameter 
Log-Normal, and Log-Pearson Type III distributions. The Three 
Parameter Log-Normal distribution was selected as providing the best 
fit. 

Flood frequency curves for each station are presented in Figures 4-1 
through 4-7. Upper and lower 95% confidence limits are shown on these 
curves. 

The flow regime on the Grand, Conestogo, and Speed Rivers is greatly 
affected at the present time by the Shand, Conestogo, and Guelph Dams. 
When the frequency of flood flows of the existing river system is 
required, it is necessary to consider the effect of the existing 
reservoir system. 

The existing reservoir operating procedure was modelled and applied in 
such a way that the effect of the reservoir system is considered in 
the consistent manner required for frequency analysis (see 1.0 
Background). This was done using HEC-5, "Simulation of Flood Control 
and Conservation Systems"; a computer program prepared by the 
Hydrologic Engineering Centre, U.S. Army Corps of Engineers. As 
applied to this study, the model simulates the use of available flood 
storage in a system of reservoirs to minimize flood flows at selected 
downstream locations. 

The initial water level assumed at the reservoirs was the April 1st 
target level set out in the existing operating policy. Storage 
available for flood control at this level is as follows: 

Bel wood - 10,000 acre-feet 
Conestogo - 8,400 acre-feet 

Guelph - 4,500 acre-feet 

The effect of this flood control storage was determined for the floods 
having 5,10,20,50, and 100-year return periods with no reservoirs in 
place. The resulting flood flows were assumed to have the same return 
periods. 

The methodology used and the assumptions made are discussed in more 
detail in Technical Report #38, "Reservoir Operations". 

For completeness, the flood frequency curves for flows regulated by 
the existing reservoir system are presented in Figures 4-8 through 
4-12 for the following stations: 

Grand River at Shand Dam 
Conestogo River at Conestogo Dam 
Grand River at Gait 
Grand River at Brantford 
Speed River below Guelph 
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Figure 4-2. CONESTOGO RIVER AT CONESTOGO DAM 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 



Natural Flow with 95% Confidence Limits 
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Figure 4-3. GRAND RIVER AT GALT 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 



Natural Flow with 95% Confidence Limits 
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Figure 4-4. GRAND RIVER AT BRANTFORD 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 



Natural Flow with 95X Confidence Limits 
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Figure 4-b. SPEED RIVER BELOW GUELPrf 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1950-1979 



Natural Flow with 95% Confidence Limits 
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Figure 4-6, nith river at new Hamburg 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1951-1979 



Natural Flow with 95% Confidence Limits 



kx) ooo 



•j'i 'i < 'J 6 



!- I F 



ProbaLihly of Exceedence (percent) 
Vir. -JO 6(3 /() (iO 50 40 30 20 



I 
i 



— , . .. } ■■■ -t~ 



._.:_. 



I 

KJOJO 



1000 



i : 



? J . ..; 



--;-! 



Z2IZ2 ] .! .' 



Iffl 



mmm 



£ 1 5 .ij Ul(li6 Oi'.l 




3 ; 



j 



Ftecurrence Interval (years) 



Figure 4-7. NITH RIVER AT CANNING 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 



Natural Flow with 95% Confidence Limits 
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Figure 4-8. GRAND RIVER A[ SHAND DAM 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 



Natural Flow 
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Figure 4-9. CONESTOGO RIVER AT CONESTOGO DAM 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 
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Figure 4-1U. UKAND RIVER AT GAL I 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 
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Figure 4-11. GRAND RIVER AT BRANTFORD 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948- 1979 
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reservoir system (HEC b 
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Figure 4-12. SPEEU K1VER BELOW GUELPH 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1950-1979 



Natural Flow 



Regulated by existing 

reservoir system (HEC 5^ 
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APPENDIX A 



MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 



-A2- 



TABLE A-l MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 





SHAND DAM INFLOW 






MAX. DAILY 


MAX. INST. 


1942 






N = 23 


1943 






x = 6080 


1944 






y = 9450 


1945 






y = 276 + 1.510x 


1946 








1947 


16900 






1948 


25800 




1949 


5030 


7870 




1950 


8030 


12400 




1951 


6580 


10200 




1952 


5460 


8520 




1953 


2950 


4730 




1954 


6380 
5710 


9900 




1955 


11700 




1956 


8410 
2790 


13300 




1957 


4490 




1958 


1990 
5760 


3280 




1959 


8040 




1960 


8540 


11200 




1961 


2440 


3460 




1962 


5010 


10400 




1963 


6680 


9680 




1964 


3140 


4630 




1965 


5100 


5980 




1966 


2240 


3220 




1967 


4590 


8510 




1968 


3720 


5580 




1969 


3320 


5710 




1970 


5610 


8300 




1971 


7270 


8290 




1972 


9930 


14100 




1973 


4550 


7130 




1974 


8580 


14800 




1975 


8410 


15600 




1976 


7850 


14100 




1977 


9010 


13100 




1978 


6610 


8120 




1979 


7310 


12500 





Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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Maximum 2 -Day Mean (IOOO cfs] 
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TABLE A-2 MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 







CONESTOGO DAM INFLOW 








MAX. DAILY MA) 
7960 


(. INST 




1948 


12600 




1949 




4050 


7900 




1950 




5320 


9410 




1951 




3390 


7110 




1952 




3080 


6750 




1953 




1300 


4630 




1954 




5720 


9890 




1955 




1580 


4960 




1956 




5520 


9650 




1957 




1290 


4620 




1958 




353 
3490 


3500 




1959 


6350 




1960 




9950 


13000 


1961 




1910 


3650 


1962 




6460 


12800 


1963 




7840 


11700 


1964 




3740 


8340 


1965 




5080 


15600 


1966 




1670 


3130 


1967 




4790 


12000 


1968 




3580 


7440 


1969 




3410 


8280 


1970 




4220 


6720 


1971 




4540 


5250 


1972 




5200 


8590 


1973 




3830 


8990 


1974 




4750 


7670 


1975 




12300 


20800 


1976 




11600 


17200 


1977 




12930 


15320 


1978 




6520 


8460 


1979 




8510 


13600 


N = 


21 






x = 


6020 






y = 


10200 






y = 


3080 + 1.1 90x 









Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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Figure A-2. 
Conestogo Dam 
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TABLE A-3 MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 
NATURAL FLOW-GRAND RIVER AT GALT 







MAX. DAILY 


1943 




28700 


1944 




17600 


1945 




25600 


1946 




20200 


1947 




37200 


1948 




44100 


1949 




23200 


1950 




27600 


1951 




23400 


1952 




16600 


1953 




10700 


1954 




18500 


1955 




10500 


1956 




24900 


1957 




9610 


1958 




5850 


1959 




27600 


1960 




38600 


1961 




8500 


1962 




23100 


1963 




25400 


1964 




11500 


1965 




20100 


1966 




7100 


1967 




28300 


1968 




13100 


1969 




17500 


1970 




13600 


1971 




20300 


1972 




27800 


1973 




22600 


1974 




26500 


1975 




31500 


1976 




30800 


1977 




32800 


1978 




24100 


1979 




23000 


N 


23 




x = 


22100 




y = 


28700 




y = 


1.317 x 


-499 



MAX. INST. 



37300 
22700 
33200 
26100 
48500 
57600 
30100 
35900 
30300 
21400 
13600 
23900 



17200 
31100 




32900 
46300 
10400 
29700 
28200 
15000 
26700 
9380 
37000 
14600 
23800 
17900 
21400 
32400 
28200 
33400 
58700 
42800 
40400 
26600 
36000 



Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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Figure A- 3. 
Grand River at Gait 

Natural Flow 
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TABLE A-4 MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 



NATURAL FLOW-GRAND RIVER AT BRANTFORD 





' 


MAX DAILY 


1948 




52900 


1949 




24600 


1950 




45700 


1951 




30400 


1952 




20700 


1953 




14400 


1954 




32300 


1955 




23300 


1956 




38500 


1957 




16000 


1958 




7800 


1959 




31000 


1960 




45400 


1961 




12500 


1962 




24700 


1963 




32600 


1964 




13000 


1965 




27700 


1966 




7320 


1967 




20700 


1968 




16910 


1969 




20000 


1970 




12700 


1971 




22400 


1972 




33300 


1973 




25400 


1974 




28500 


1975 




49600 


1976 




44400 


1977 




50000 


1978 




30100 


1979 




47100 


N 


23 




X 


28600 




y = 


34200 




y » 


2670 + 


1 . 104x 



MAX. INST. 



61100 
29800 
53100 
36200 
25500 
18600 
38300 



26000 
40500 




35700 
55500 
14100 
34600 
36900 
17000 
33600 
9320 
36500 
18400 
24300 
17500 
24300 
38000 
33800 
36000 
64000 
52200 
53200 
33100 
52700 



Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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Figure A-4. 
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TABLE A-5 MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 



NATURAL FLOW-SPEED RIVER BELOW GUELPH 





MAX. DAILY 


1950 


_....— 


1951 


3000 


1952 


1580 


1953 


1390 


1954 


2460 


1955 


2060 


1956 


3260 


1957 


1210 


1958 


550 


1959 


2940 



N = 19 

x = 2460 

y = 2870 

y = 91.5 + 1.128x 

MAX. INST. 

6720 
3780 




1960 

1961 

1962 
1963 
1964 
1965 
1966 
1967 
1968 
1969 
1970 
1971 
1972 
1973 
1974 
1975 
1976 
1977 
1978 
1979 



2860 

1190 

1980 
2840 

870 
2760 

820 
3780 
1330 
2520 
1400 
1530 
3480 
2670 
3310 
3370 

3400 
2010 
2800 



3010 



1910 



2170 
3000 
1020 
3460 
920 
4810 
2300 
2750 
1640 
1700 
3730 
2930 
3430 
4480 
3500E 
3700 
2280 
3340 



Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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Figure A-5. 

Speed River at Guelph 

Natural Flow 
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Maximum Daily Mean (1000 cfs) 
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TABLE A-6 MAXIUMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 



NITH RIVER AT NEW HAMBURG 







MAX. DAILY 


1951 




3490 


1952 




2750 


1953 




2620 


1954 




4470 


1955 




3250 


1956 




4990 


1957 




2920 


1958 




2300 


1959 




5090 


I960 




7130 


1961 




2180 


1962 




5480 


1963 




7520 


1964 




2700 


1965 




5680 


1966 




1640 


1967 




7460 


1968 




5140 


1969 




5080 


1970 




3780 


1971 




3030 


1972 




5090 


1973 




4680 


1974 




6920 


1975 




10700 


1976 




9380 


1977 




9110 


1978 




4910 


1979 




9390 


N 


22 




x = 


5860 




y = 


8260 




y = 


1:429 x ■ 


• 113 



MAX. INST. 



5570 
3630 




7650 




5990 
8820 
2470 
7010 
9220 



3750 



9450 



2230 



11800 

9880 

6950 

4440 

3680 

6290 

6800 

11000 

16000 

11400 

12800 

7030 

13700 



Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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TABLE A-7 MAXIMUM INSTANTANEOUS FLOWS VERSUS DAILY MEAN FLOWS 



NITH RIVER AT CANNING 









MAX. DAILY 




1948 




10500 




1949 




6130 




1950 




11600 




1951 




5700 




1952 




3900 




1953 




2680 




1954 




10600 




1955 




7110 




1956 




8100 




1957 




4050 




1958 




1980 




1959 




6000 




1960 




8760 




1961 




2920 




1962 




5910 




1963 




8590 




1964 




2560 




1965 




8380 




1966 




2500 




1967 




7440 




1968 




5500 




1969 




5860 




1970 




4550 




1971 




4070 




1972 




6600 




1973 




5260 




1974 




8460 




1975 




11300 




1976 




9770 




1977 




10900 




1978 




6180 




1979 




10900 


N 


= 


29 




X 


= 


7040 




y 


= 


8830 




y 


= 


634 + 


1.164x 



MAX. INST 

14900 
9910 

12600 
6720 
3970 
4590 

12900 
8940 
9420 
4850 
2360 
6510 

10200 
6500 
7320 

10300 



3610 



12500 



3540 



10400 



7040 



7010 

5100 

4280 

7350 

7160 

9710 

14800 

11100 

12800 

7900 

13500 



Maximum instantaneous flows enclosed in boxes are estimated 
based on corresponding maximum daily flows. 
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APPENDIX B 



MAXIMUM NATURAL FLOW VERSUS CORRESPONDING OBSERVED FLOW 
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TABLE B-l 


MAXIMUM 


NATURAL FLOW VERSUS CORRESPONDING 








OBSERVED FLOW 










GRAND RIVER AT GALT 


CORRESPONDING 








MAXIMUM NATURAL 


OBSERVED FLOW 


YEAR 


DATE 
Apr 2 


FLOW PEAK (cfs) 
37300 


PEAK (cfs) 


1943 




1944 


Mar 


25 


22700 


28900 


1945 


Mar 


16 


33200 


19200 


1946 


Mar 


7 


26100 


30200 


1947 


Apr 


12 


48500 


41400 


1948 


Mar 


20 


57600 


46300 


1949 


Mar 


23 


30100 


30400 


1950 


Apr 


5 


35900 


38600 


1951 


Mar 


31 


30300 


23800 


1952 


Apr 


2 


21400 


17300 


1953 


May 


3 


13600 


15000 * 


1954 


Apr 


17 


23900 


20200 * 


1955 


Mar 


12 


17200 


15400 


1956 


Apr 


5 


31100 


30800 


1957 


Feb 


27 


12200 


14900 * 


1958 


Mar 


29 


7210 


4190 


1959 


Apr 


6 


32900 


20100 


1960 


Apr 


4 


46300 


28600 


1961 


Feb 


25 


10400 


6680 


1962 


Mar 


30 


29700 


13800 


1963 


Mar 


28 


28200 


14400 


1964 


Apr 


8 


15000 


5730 


1965 


Apr 


3 


26700 


22100 


1966 


Mar 


19 


9380 


4590 


1967 


Apr 


3 


37000 


30800 


1968 


Mar 


21 


14600 


7010 


1969 


Apr 


19 


23800 


24400 


1970 


Apr 


15 


17900 


8830 


1971 


Apr 


13 


21400 


13300 


1972 


Apr 


19 


32400 


30100 


1973 


Mar 


12 


28200 


17300 


1974 


May 


17 


55200 


54600 


1975 


Apr 


19 


58700 


30100 


1976 


Mar 


21 


42800 


20900 


1977 


Mar 


13 


40400 


20100 


1978 


Apr 


12 


26600 


16400 


1979 


Apr 


14 


36000 


27900 



Flows printed in bold are estimated based on daily mean flows 
* Flow is not maximum observed in this year. 
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TABLE B-2 MAXIMUM NATURAL FLOW VERSUS CORRESPONDING 









OBSERVED FLOW 








GRAND 


RIVER AT BRANTFORD 
MAXIMUM NATURAL 


CORRESPONDING 
OBSERVED FLOW 


YEAR 


DATE 


FLOW PEAK (cfs) 


PEAK (cfs) 


1948 


Mar 


20 


61100 




1949 


Mar 


23 


29800 


40000 


1950 


Apr 


5 


53100 


25700 


1951 


Mar 


31 


36200 


12900 


1952 


Apr 


2 


25500 


20800 


1953 


May 


3 


18600 


17200 * 


1954 


Apr 


17 


38300 


39000 * 


1955 


Mar 


12 


26000 


25800 


1956 


Apr 


5 


40500 


40300 


1957 


Feb 


27 


20300 


20600 * 


1958 


Mar 


29 


11300 


6270 


1959 


Apr 


6 


35700 


25800 


1960 


Apr 


4 


55500 


40600 


1961 


Feb 


25 


14100 


14600 


1962 


Mar 


30 


34600 


19100 


1963 


Mar 


28 


36900 


24300 


1964 


Apr 


8 


17000 


7900 * 


1965 


Apr 


8 


33600 


29000 * 


1966 


Mar 


19 


9320 


2960 


1967 


Apr 


3 


36500 


34900 


1968 


Mar 


21 


18400 


12700 * 


1969 


Apr 


19 


24300 


25800 


1970 


Apr 


15 


17500 


12800 


1971 


Apr 


13 


24300 


16500 


1972 


Apr 


19 


38000 


35700 


1973 


Mar 


12 


33800 


23700 


1974 


May 


17 


64700 


63000 


1975 


Apr 


19 


64000 


37200 


1976 


Mar 


21 


52200 


28200 


1977 


Mar 


13 


53200 


32400 


1978 


Apr 


12 


33100 


21100 


1979 


Apr 


14 


52700 


46400 



Flows printed in bold are estimated based on daily mean flows 
* Flow is not maximum observed in this year. 
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FIGURE B-3 


MAXIMUM 


NATURAL FLOW VERSUS CORRESPONDING 








OBSERVED FLOW 










SPEED RIVER 


1 BELOW GUELPH 


CORRESPONDING 








MAXIMUM NATURAL 


OBSERVED FLOW 


YEAR 


DATE 


FLOW 


PEAK (cfs) 


PEAK (cfs) 


1948 












1949 


Mar 


23 








1950 


Apr 


5 




6720 




1951 


Mar 


31 




3780 




1952 


Apr 


2 




1870 




1953 


May 3 




1660 




1954 


Apr 


V 




2870 




1955 


Mar 


12 




2420 




1956 


Apr 


5 




3730 




1957 


Feb 


27 




1460 




1958 


Mar 


29 




712 




1959 


Apr 


6 




3410 




1960 


Apr 


4 




3010 




1961 


Feb 


25 




1430 




1962 


Mar 


30 




2170 




1963 


Mar 


28 




3000 




1964 


Apr 


8 




1020 


. 


1965 


Apr 


8 




3460 




1966 


Mar 


19 




920 




1967 


Apr 


3 




4810 




1968 


Mar 


21 




2300 




1969 


Apr 


19 




2750 




1970 


Apr 


15 




1640 




1971 


Apr 


13 




1700 




1972 


Apr 


19 




3730 




1973 


Mar 


12 




2930 




1974 


May 


17 




5850 




1975 


Apr 


19 




4480 




1976 


Mar 


21 




3500 


1770 


1977 


Mar 


13 




3700 


1830 


1978 


Apr 


12 




2280 


1220 


1979 


Apr 


14 




3340 


2020 



Flows printed in bold are estimated based on daily mean flows 
* Flow is not maximum observed in this year. 
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Appendix C 

FREQUENCY DISTRIBUTIONS 



Excerpted from "Flood Damage Reduction Program Flood Frequency Analysis", 
R. Condie, G.A. Nix, and L.G. Boone, Environment Canada, January 1977. 



This appendix gives a brief mathematical summary of the frequency 
distributions used in the program. Only those features of the distributions 
that are necessary for an understanding of the program are given, along 
with the solution. Four frequency distributions are considered, the Gumbel 
I, the Log-Normal, the Three Parameter Log-Normal and the Log-Pearson III, 
and as far as possible they are fitted to the flood data by the method of 
maximum likel ihood. 

For the Three Parameter Log-Normal with the probability density function in 
its conventional form, there will seldom be a maximum likelihood fit when 
the untransformed data exhibit negative skewness and in such cases there 
can never be a moment fit either. On the other hand, a maximum likelihood 
fit may not be obtained when the data have a very small positive skewness 
but a moment fit fs always possible. In those cases, the program switches 
to a moment fit and identifies it as such on the output. 

The Log-Pearson III will be bounded above when the logarithms of the flood 
data h3ve neaativc- skewness; wnen positive the distribution will be bounded 
below. In general, for Canadian rivers, this distribution will be bounded 
above and it is difficult to interpret "the analysis. Sometimes the maximum 
likelihood fit gives results which seem, subjectively, better than those 
obtained by iroments, and vice versa. For this distribution the frequency 
regime is given for both maximum likelihood and moment fitting. 

GENERAL: 

Population Statistics 

Given a sample Xj, Xj>, ... xN, the best estimates of the population mean, 
standard deviation, and skew and kurtosis coefficients as computed herein 
are given by 

Mean, x * (l/N)Ex G(l) 

Standard Deviation, 5 * ((l/(N-i)]I(x-x) 2 } G(2) 

Skew Coefficient , t l m (N 7 / [ (N-l) (N 2) ] } (> 3 /S 3 ) C(3) 

Coefficient of Kurtosis. Z z * ( [N 2 CN*1 ) ]/ [ (N- 1 ) (N- Z) (N-3) ] } (m^/5 1 *) G(4) 

where m„ and m. , the third and fourth central ir.orncr.es arc defined by 
3 *• 

m 3 «= Cl/N)£(x-x) 3 G(5) 

and m. ■ (l/K)£(x-x)* C(6) 
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The summations are carried out over the N terns of the data series. When 
working in logarithmic units, the same definitions apply with x replaced by 
its natural logarithm, lnx. There is no universally accepted system of 
statistical notation and for that reason symbols are identified throughout 
the text. Distribution parameters when estimated from a sample are 
distinguished by putting a circumflex over them. 

The Cumultive Probability Function or General Frequency Equation 

This function gives the probability of an event less than a stated 
magnitude. For most frequency distributions it cannot be expressed in 
closed form. In this program which deals specifically with flood frequency 
analysis it is rearranged to give the T year flood as a function of the 
distribution parameters only, or a function of these parameters and the 
standard normal deviate t. 



Maximum Likelihood Estimators 

or a given probability density function * (x; a, b, c, ...), where a b 
, etc. are parameters, the joint probability that a sample of N values xi , 
2,...xN came from that distribution is l 

L = W*tJ a » *>. c, ...)4»(x 2 ; *c f b f c, ...) ... $(x N ; a, b, c, ...) 



For 

c 

x 2 



L is the likelihood and the principle of maximum likelihood states that the 
values of a, b, c, etc. should be chosen to maximize L. When dealinq with 

WilSJ*-? 1 ^ ?£'?"! 13 i SUdlly hecomes simpler t0 inax "nize the natural 
logarithm of the likelihood function InL. 

The system of equations which when solved achieve this maximization are the 
Maximum Likelihood Estimators for the parameters of the distribution, and 
tne system is 

91nL/Sa ■ 31nL/3b * 31nL/3c etc. = 
Minimum variance estimates of parameters are obtained by this method. 

The Asymptotic Standard Error of Estimate 

If z is a function of sc»y three statistical variables a, b and c, each 
subject to sampling variance, then the variance of z is 

var(z) * (3z/3a) 2 var(a) + C3z/3b) 2 var(b) * (3z/3c) ? var(c) 

4 2(3z/3i)(Dz/3b)cov(a,b) * 2(3z/3a) (3z/3c)cov(a,c) 

+ 2(3z/3b)(3z/Dc)cov(b,c) C(5) 
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Equation G(5) is the general variance of estimate equation and is 
applicable to a function of any number of variables. The required oartial 
derivatives are obtained from the general frequency equation and a^e then 
evaluated at the estimated values a\ d*, t. The required variances and 
covariances are obtained from the inverse dispersion matrix using the 
estimated values -a\ 'fr, t". Sometimes, particularly fcr two parameter 
distributions Equation 3(5) reduces to a relatively simple function and the 
Asymptotic Standard Error of Estimate is root var(z). For the three 
parameter distributions the expression is usually too complex to be shown 
within the confines of a single page. In this text for the Three Parameter 
tog-Normal and Log-Pearson Type III distributions only the required oartial 
derivatives and the inverse dispersion matrix are given. Matrix invasion 
to give the required variances and covariances is performed by a subroutine 
in the program. 

THE DISTRIBUTIONS USED 

A. THE GUMBEL I DISTRIBUTION: 

Fisher and Tippett's first asymptotic distribution of extreme values, 
sometimes known as the double expotential distribution, was extensively 
investigated by E.J. Gumbel and proposed for use in flood frequency 
analysis. In view of his many contributions in* this field it is usually 
known in hydrology as the Gumbel I distribution. 

It has two parameters, scale and location, the va^iate is unbounded either 
above or below, and has fixed coefficients of skew and kurtosls of 1.14 and 
5.4 respectively. Because of its simplicity and also having some 
theoretical justification, use of this distribution was formerly much 
favoured in flood frequency analysis. The lack of a lower boundary can at 
times lead to negative floods at shorter return periods, a condition which 
is not physically possible. On the assumption that the probability of 
occurrence of less than zero is the probability mass of zero flood, the 
implicaton is that the river may be dry for an entire year, if it is the 
annual flood series which is being analysed. 

8y considering these few simple characteristics, a reasonable assessment 
can be made of how well this distribution has been fitted to the flood data 
being analysed. The occurrence in the analysis of negative values at 
shorter return periods should be examined jointly with the plotten oata and 
of the population coefficients of skewness and kurtosis as computed from 
the sample, and if their values are fairly near the theoretical values of 
1.14 and 5.4 it can reasonably be concluded that a good fit nas been 
obtained. 

The Probability Density function 



, . -a(x-u) 
♦ (x) - .e- a ' X - u V e A(l) 
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'he Cumulative Probability Function 
-a(x-u) 
P - e A(2) 

where P is the probability of a value less than x, a is a scale parameter, 
a > o, and u is a location parameter. These parameters are estimated from 
the sample using maximum likelihood theory. To distinguish these estimated 
parameters they eire designated a and u. 

Maximum Likelihood Estimators 

Applying maximum likelihood theory to the p.d.f. at Equation A(l) qives the 
following * ' ■ 

ainL/aa ■ N/a - E(x-u) * Z{(x-u)exp. - a(x-u) > » A(3) 

31nL/3u ■ Na - a£{exp. - a(x-u)} * X(4) 

where the summations ire taken over the N terms of the sample. 
Re-arranging Equations A(3) and A(4) gives. 

a£xe" ax - (ax-l)2e~ ax - A(5) 

u ■ (1/aKlnN-lnEe" *) A(6) 

Equations A(5) and A(6) were given by Kimball, Reference (5). The curve 
corresponding to Equation A(5) has two inflexions and is solved by using 
the Bolzano method first to isolate the root clear of the inflexions, then 
by the Newton-Raphson method. The estimated parameter a is then 
substituted in Equation A(6) to give ft. 

The Frequency Regime 

Re-arranging Equation A(2) and replacing a and u by their estimates. 3. 0". 
gives ■ • 

x t « Q - (l/a)ln.ln{T/(T-l)} A(7) 

where x t is the value of the variate corresponding to a return period of T 
years. ^ 

The Asymptotic Standard Error of Estimate 

The inverse dispersion matrix can be derived from Equations A(3) and A(4) 
Inverting, simplifying, and substituting in the general equation for the 
variance of estimate gives 

var(x t ) - {l/(a 2 N)Hl*(l-Y+y t ) 2 /(n 2 /6)> A(8) 

where r is Euler's constant, approximately 0.5772 and y* is the reduced 
variate a(x t -u). 

This expression is given by E.J. Gunbel, Reference (6) following the work 

jjf B.F Kimball, Reference (7). Replacing a and u by their estimates 2 and 

-j, and taking the square root of var(x t ) gives the Asymptotic Standard 
Er>-or of the T year event. 
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B. THE LOG-NORMAL DISTRIBUTION: 

Hydrologic events can seldom be described by the normal distribution since 
these events are more commonly skewed.. Many theoretical transformations 
have been developed for normalizing skewed distributions, because it is 
generally easier to draw statistical conclusions from the normal 
distribution whose exact theoretical properties are well known. In flood 
frequency analysis it is generally necessary to extrapolate far bevond the 
range of observed events and the idea of a normalizing transformation is 
very attractive. One such transformation leads to the Log-Normal 
distribution. This distribution has been around for quite a long time. 
Galton noted in 1875 that the geometric mean of some types of data series 
was a better descriptor than the arithmetic mean. The geometric mean is of 
course the arithmetic mean of the logarithms, or rather the ant i log of 
such, and hence the name Log-Normal. Some older text books still refer to 
this distribution as the Galton distribution. Hazen in 1930 applied this 
distribution to flood frequency analysis and the distribution is known by 
three names; the Log-Normal, Hazen, or Galton distributions. Essentially 
use of this distribution implies that the logarithms of the basic data to 
be analysis are normally distributed. 

It is a simple two parameter distribution and the two parameters jointly 
influence the location of the mean, the scale, and the shape. The v.'.riate 
is bounded below at zero and has no upper boundary, and this distribution 
cannot have negative skewnass. This latter feature is very important 
because the Log-Normal distribution can he fitted to negatively skewed data 
but the analysis will be grossly in error. By virtue of the logat'».n,..ic 
transform the theoretical coefficients of skew and kurtosis of the 
transformed data are respectively 0.0 and 3.0. 

These theoretical characteristics give a handy subjective assessment of the 
validity of fitting the Log-Normal distribution to the data sample of 
annual floods. The program output gives best estimates of the coefficients 
of skew and kurtosis of the logarithms of the population as estimated from 
the available sample of annual floods. If their values do not vary wry 
much from the theoretical values of 0.0 and 3.0, then the distribution will 
give a fairly good estimate of the flood frequency regime at the 
particular location. If the skewness of the logarithms is substantially 
positive then floods at the longer return periods will be cnderestuiated 
and if negative, tnese floods will be over estimated. These conditions can 
be verified by examining the plotted data. 

Experience with Canadian rivers shows that the logarithmic transform of the 
annual flood series generally results in overdoing the normalizing 
process. Basic data which are positively skewed show sometimes a auite 
substantially negative skew coefficient of the logarithmic transform. 
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The Probability Density Function 

xo^T 2o- B(l) 

where m and o 2 are respectively location and scale parameters of the 
transformed variate lnx. For the untransformed variate they jointly 
influence the skewress, but the lower boundary is fixed at zero. 

The Cumulative Probability Function 

The cumulative probability function cannot be expressed in closed form and 
instead, the general frequency equation is used 

H s exp. (nt+to) B(2) 

The parameters are replaced by their estimates m and $ and I is the 
standard normal deviate at. the required return period. 

Maximum Likelihood Estimators 

Applying maximum likelihood theory to the p.d.f. at Equation 8(1) gives 

»lnL/3a = -(Nn)/cx 2 ♦ l/a 2 Zlnx « " B(3) 

UnL/3(a 2 ) - -N/(2o 2 ) * l/(2o ,| )I(lnx-») 2 - B(4) 

Re-arranging Equations B(3) and B(4) gives 

» » (l/N)Hnx B(5) 

a 2 « (l/N)E(lnx-n0 2 B(6) 

Equation B(6) gives a biassed estimate of o 2 and is corrected for the bias 
to give 

o 2 • (l/(N-n}J:nnx-m) 2 B(7) 

Equation B(5) and B(7) give minimum variance unbiassed estimtes ft and o 2 
and m and o. 

The Asymptotic Standard Error of Estimate 

The inverse dispersion matrix can be derived from Equations 8(3) and 3(4). 
Inverting, simplifying, and substituting in the general equation for the 
variance of estimate gives 

var(x t ) - x t 2 (a 2 /N)Cl** 2 /2) B(8) 

Replacing o by its estimte § and taking the square root give the Asymptotic 
Standard Error of Estimate of the T yetr event. 
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C. THE THREE PARAMETER LOG-NORMAL DISTRIBUTION: 

The success of any normalizing transformation can be measured by computing 
the coefficients cf skew and kurtosis of the transformed data, and tne 
closer these coefficients are to 0.0 and 3.0, the more successful the 
transformation is. The straightforward transformation used in the 
Log-Normal distribution is not markedly successful in this respect and it 
has been found that the inclusion of a third parameter will greatly improve 
the transformation, hence the name, the Three Parameter Log-Normal 
distribution. In its most simple terms, the function y = ln(x-a) is 
distributed normal ly. 

The three parameters of the distribution jointly influence the location of 
the mean, the scale and the shape. The variate is bounded below by the 
parameter a, and has no upper boundary and in its conventional form the 
distribution cannot have negative skewness. When developing this program 
some theoretical work was done on the negatiely skewed Three Parameter 
Log-Normal distribution. The variate in that distribution is unbounded 
below and has a positive upper boundary at the parameter a. Although this 
form cf the distribution was successfully programmed, it was generally 
found that the upper boundary parameter was very close to the maximum 
observed flood in the data series. Subjectively, this form of the 
distribution should not be used for flood frequency analysis. Further 
development work compared the results obtained by fitting the distribution 
by moments with those obtained by a maximum likelihood fit. The maximum 
likelihood fit was clearly superior in that it gave a smaller standard 
error of estimata of the T-year event, and also the computed coefficient of 
skew of the transformed data was much closer to the theoretical value of 
0.0. For these reasons, only the maximum likelihood fit is used, except in 
the circumstances mentioned earlier, when a moment fit is used. 

The Three Parameter Log-Normal is an exceptionally flexiole distribution 
and \/ery well suited to flood frequency analysis. When the parameter a 
becomes zero, then the distribution becomes the straightforward Log- Normal , 
and it can also be shown theoretically that tne Gumbel I distribution is a 
special case of the Three-Parameter Log-Normal. Extensive testing during 
the development of this program showed that the coefficient of skew of the 
transformed data averaged 0.055, with a low value of 0.0008, an excellent 
check on the success of the fit. 

Like the Log-Normal, the validity of fitting the Three Parameter 
Log-Normal distribution can easly be assessed by comparing the computed 
values of the skew and kurtosis coefficients with their theoretical 
values. One minor drawback is that where the Log-Normal distribution shows 
skewness of the transformed data very close to 0.0, the Three Parameter 
Log-Normal will usually improve on this but at the expense of an increased 
standard error ef estimate. This is caused by the introduction of the 
third parameter which is itself subject to sampling variance and is 
reflected in the overall standard error of estimate. This is not ^ery 
serious, however, since the T year event obtained from either of these 
distributions in such cases is pretty well the same. 
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The Probability Oensity Function 

♦ (x) • L— exp. - Hn(x-a) - m) 2 

where m and o^ are respectively location and scale parameters for the 
transformed variate In(x-a). For the un trans formed variate they jointly 
influence the skewness and a is the lower boundary of the variate x. 

The Cumulative Probability Fucntion 

Like the log-normal distribution, the cumulative probability function 
cannot be expressed in closed form and the general frequency equation is 

x t • a ♦ cxp(m+to) C(2) 

The parameters are replaced by their estimates a\ m* and 6 s and t is the 
standard normal deviate at the required return period. 

Maximum Likelihood Estimators 

Applying maximum likelihood theory to the p.d.f. at Equation C(i) gives 

81nL/8a » (l-m/a 2 )Z(x-a)- 1 + (l/o 2 )E{[ln(x-a) ](x-a)" 1 } » C(2) 

ainL/fcn « (l/o 2 )E{ln(x-a)-m} « C(3) 

*lnL/3(o 2 ) = -n/(2o 2 ) ♦ l/(2o'»)E{lnCx-a)-a} 2 - C(4) 

Equations C(2), and C(3) and C(4) can be re-arranged to give an equation in 
a only f(a) = 

E(x-a)- 1 {(l/N)rin 2 (x-a) - [(l/N)Zln(x-a) ] 2 - (?/N)Eln(x-a) } 

+ E{ln(x-a)>(x-a)- 1 » C(5) 

When developing a solution for Equation C(5) it was noticed that f(a) and f 
(a) were negative near x min., with both increasing as a decreased. A 
Newton-Raphson solution was used with iteration starting at O.Sx min. and x 
min. To avoid missing such a root, f(a) is evaluated at these points and 
checked for a sign change, and a Solzano method is included in the program 
to isolate any such roots. 

Having found a, then 

m « (l/N)IlnCx-a) C(6) 

and, correcting for bias 

o 2 - {l/(N-l)}L{lr.(x-a)-m} 2 (7 ) 

Hence the estimated parameters a, <ti and $2. 
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The Asymptotic Standard Error of Estimate 

For the general variance of estimate equation, the required partial 
aenvatives are 

'V Ja - ' C (8) 



<s 



C(9) 
C(10) 



C(ll) 



-___ Jx t /8n • exp. (ra+ta) 

— *- >x t /3ra 2 = (l/2o 2 )exp.(m*ta) 
and the inverse dispersion matrix V~* is 

(l*o 2 )exp. (2o 2 -2m) exp. (o 2 /2-m) -exp. (o 2 /2-m) 

N/o 2 exp.(o 2 /2-m) i Q 

-exp. (o 2 /2-n) i/ 2ff 2 

Equations C(8) through C(10) are evaluated at the estimated values m and 3 
with the appropriate t for the required probability level. The elements of 
the matrix are evaluated at tf and fi 2 and inverted-to give the variances and 
coyanances. Substitution in the general variance of estimate eouaticn and 
the 1 ? 9 ear eJent* r ° 0t ^ ^ Aswtotic Stand ^d Error of Estimate of 

Fitting by Moments 

Taking moments of the p.d.f. given by equation C(l), re-arranging and 
replacing them by their sample estimates gives 

k 3 ♦ 3k - gj = 

A - ln(SA) - l/21n(k 2 +l) C(14) 

o 2 - ln(k 2 +l) C(15 , 

«nere s and gj are computed from Equations G(2) and G(3). 

For positive skew and from Descartes' rule o* signs Equation C'12) has only 

4V!iiJ *°ni^ \1* positive ' The cubic ** solved directly for k, then 
a", to and o -ollow. For neaative skew, Equation C(12) has only one real 
root, which is negative. Equation C(14) therefore cannot be evaluated. 

The Asymptotic Standard Error for Estimate of the T year event is computed 
according to Kite, Reference (8). The method is quoted here practically in 
full. Parameterization is slightly different to avoid confusion f'ie T 
?Sf r »H en i C *!) ?! ^Pressed in terms of the first moment about the onoin 
the second and third central moments of the population and the standard ' 

ske r w a coe?f!c'e S nt"by Den ? tln9 ^ """^ ** "^ ^ ^ ^ the W*l«"«" 

1/2 
\ * » " VI ♦ Fy 2 C(16) 

where F ■ cxp. ([lnf Hk-) ] l 2 t - [ Inf \*k 7 ) ]/2} - 1-0 C(171 

k 



C(12) 
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and k is the root of Equation C(12). 
The direct solution of Equation C(12) is 

k « (l-w 2 /3)/ w l/3 

C(18) 

where w * ((yi ♦4) 1 / 2 - v, }/2 

C(19) 
Since Xl = "3 /u 2 3 / » Equation C(16) can be rewritten 

x t = z = f(u lt u 2 , u 3 , t ) 
t is invariable and from Equation G(5) 

var(z) = Oz/a^jivarC^) + <3z/3y 2 )2 V ar(y 2 ) * (3 2 /3M 3 )var(y 3 ) 

♦ ZCSz/SM'pcaz/a^covCu'^up + 2(3z/3y' i )C3z/3 l : 3 )cov(y',y 3 ) 

♦ 2(3z/3y 2 )Oz/3y 3 )cov(y 2 , U3 ) C(2Q) 

The required partial derivatives ire 

Jz/ay'j - 1 C(21) 

8z/3w 2 » (l^Uj 1 2 )(F-3 Yl 3F/3 Yl ) C(22) 

3z/3y 3 » (l/y 2 )(3F/3T 1 ) C(23) 

F has been defined at Equation C(17) and 

8F/3 Yl « (3F/3k)(3k/3w)(3w/3 Yl ) C(24) 

with 

3F/3k ll^ 2 " l|[ exp " {Iln(1+k2)ltl/2 " ( ln ( 1+k2 )3/ 2 >J # 

[a[in(i^), 1/a "7-3?] + ^* s| C(25) 

**•-*• fe^J c(27) 
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The computed values of w and k, and the appropriate t for the required 
probability level, and replacing Tfi by its sample estimate $, will then 
substituted in Equations C(25) through .C(27) give SP/W from Equation C 
(24). Substituting for F and Why*, and replacing u? ^y its sample 
estimated \ n Equations C(21) through C(23) give the required partial 
derivatives from Equation C(20). 

The required variances and covariances are 
var v " v 2 /N 



var y 2 - | (^ - »*] 



C(28) 

C(29) 

v«r u 3 - £ (w 6 - w 3 2 - 6 V2 •■ 9y^ C(30) 

cov (y'j, up * p f /|, C(31) 

cov (v'j, u 3 ) * i (^ - 3u*) C(32) 

cov ( V m 3 ) - i ( y 5 - A V ^ 2) C(33) 

For the Three Parameter Log-Normal distribution Uie population central 
moments U2 through ug are defined by 

V 2 - c° 2+2 (e a2 - 1) C(34) 

W* 4 Ce' 2 -!)* (.*•*♦ 2.*^ « 3.** B CCSS) 

. 5o 2 /2+S f 10a 2 . 6o 2 A lrt 3o 2 o 2 • 

5 te - 5e + 10e ° . I0e° ♦ 4) C(37 ) 

, -3o 2 *6 fA 15o 2 , IO0 2 . ,_ 6a 2 ,„ 3a 2 o 2 
M 6 e < e - 6e ♦ lSe - 20e ♦ is c . 5) C(38) 

Replacing o* 2 by its sample estimate t 2 in Equations C(34) through C(33) 
gives the central moments p? through p 6 and hence the variances and 
covariances from Equations C(28) through C(33). Substitution in Equation 
C(20) gives the variance of the T year event. Finally taking the square 
root gives the Asymptotic Standard Error. Occasionally this variance may 
be negative and the Asymptotic Standard Error is printed out as zeroes. 

D. THE LOG-PEARSON TYPE III DISTRIBUTION 

A group of frequency distributions can be derived from a generalized 
differential equation proposed by Karl Pearson. This generalized equation 
has four constants and by equating some to zero or to each other and 
solving the differential equation a series of syrr.etrical or skewed 
distributions are found. One solutions leads to a two-parameter 
distribution and by combinations of these parameters a variety of U and 
J-shaped distribut'ons are obtainable. Reverting to tne generalized 
differential equation, a possible solution leacs to °earson's Type III 
distribution and if it is assumed that the logarithm of the van ate follows 
a Pearson Type III distribution, the distribution of the var i ate itself is 
the Log-Pearson Type III. 
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[n 1967 the U.S. Water Resources Council recommended Us use as the 
standard method of flood frequency analysis with provision for departure 
when justified. The Council's decision was based on a need for 
standardization, convenience, since many agencies used it, and on the 
analysis of flood data from 10 American rivers. In developing this program 
and using data from 37 long term hydrometric stations on Canadian rivers it 
is unlikely that the same recommendation would be made. 

The Log-Pearson Type III distribution has three parameters and Bobee, 
Reference (1), shows the many different shapes the distribution can take, 
depending on the relationship between the parameters and their signs. Not 
all these shapes seem acceptable in flood frequency analysis and the user 
of this program is advised to study Reference (4), particularly with regard 
to the distribution shape. The range of the variate also depends on the 
sign of one of the parameters. The variate may hav,? a positive lower 
boundary and be unbounded above, or, it may have a zero lower boundary with 
a positive upper bound. A positive upper bound was found in about 67 
percent of the test stations and, depending on the position of this upper 
bound relative to the largest observed flood in the data sample, it can be 
^ery difficult lo judge how well this distribution fits the data. Cases 
have been found where an upper boundary parameter was found less than the 
largest observed flood. 

The validity of the fit of the Log-Pearson Type III distribution to the 
flood series cannot be assessed so easily as for the Gumbel I, the 
Log-Normal or the Three Parameter Log-Normal distributions although some 
assistance is given to '^ne user from the proqram output. When combinations 
of parameters lead to apparently impractical shapes of the distribution, or 
where an upper bound is found to be less than the greatest observed flood 
uhe fiood frequency regime is given in the output as a set of zeros. The 
distribution parameters will always be given in the output and the user may 
relate them to the combinations shown in Reference (4). It is vary diffi- 
cult to give any advice en the interpretation of an upper boundary and its 
implications. Points to be considered are the proximity of the upper 
Doundary to the largest observed flood and the general appearance of the 
plot. Sometimes fitting by the method of moments appears better than the 
maximum likelihood fit and for that reason the flood frequency reqimp is 
given using both fitting methods. 

The Probability Density Function 

^( X ) . exp{-flnx-m)/a) {(lnx-m)/a}b-l 
|a|xrCb) 

where a, b and n are respectively scale, shape and location parameters and 
]_ is the gamma ^unction of the argument within parenthesis. 

Bobe'e, Reference (*), shows the many different shapes which this p.d.f. can 
take, depending on the relationship between the parameters, and not. all 
these shapes are credible in flood frequency hydrology. This is particu- 
larly so in cases where "a" is negative, and th'is is Generally so for flood 
series on Canadian rivers. 



D(l) 



-C14- 



The Cumulative Probability Function 

Again, the cumulative probability function cannot be expressed in closed 
form. Making the substitution y = (lnx-m)/a in the p.d.f. at 0(1) and 
using the Wilson-Hifferty approximation for chi-square. Reference (9). 

lnx t - n ♦ a{t/(3b l/s ) - l/(9b 2/3 ) ♦ b 1 ' 3 } 3 D(2) 

where Inx^ is the logarithm of the T year event and t is the standard j 

normal deviate for the required return period. Then x^ s exp(lnx t ). 

Maximum Likelihood Estimators ! 

i 

Applying maximum likelihood theory to the p.d.f. at Equation 0(3.) gives 

. 
31nL/3a » (l/a 2 )I(lnx-m) - (l/a)Nb « 0(3) 

31nL/3b - -N*(b) ♦ Iln{(lnx-m)/a} - 0(4) 

31nL/3m » 1/a - (b- 1)1 (lnx-m) - 0(5) 

where y (b) is the Digamma or Psi function, 31nl*(b)/3b. 

Equations (3), (4), and (5) are three simultaneous transcendental eauations 

in a, b, and m, but Equation (41 can be r educed to a singla transcendental 

equation in m by re-arranging Equations (3) and (5) to give ! 

a « (l/Nb)£ln(x-m) D(C>) 

b - Klnx-mj-VCtClnx-m)- 1 - N 2 /I(lnx--n) } 0(7) 

Substitution of Equation (6) in Equation (4) eliminates a, then 
substitution of Equation (7) in Equation (4) eliminates b. The Psi 
function can be replaced by a combined recurrence formulate and asymptotic 
expansion where (b) tends to 

ln(b+2) - 1 - 1 * 1 , - 1 , - 1 -1 0(8) 

2(b+2) 12(b+2) 2 l20(b+2) 4 252(b+2)° (1+b) F 

It is quite clearly impossible to show these substitutions within the 
confines of one page. The distribution at Equation (1) can be either 
positively or negatively skewed. It can be shown quite easily by 
considering the moments of the distribution that if the skew coefficient of 
the logarithms is positive, a must be positive and therefore from Equation 
(4), (lnx-m) must always be positive for all values of Inx . It follows 
then that the distribution is boundea below ai exp. m. It is more common 
in flood series from Canadian rivers to find that the skew coefficient of 
the logarithms is negative, and (lnx-m) must always be negative. Thus the 
distribution is bounded above at exp. m. 
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The solution of Equation (4) with the necessary substitutions can therefore 

thP \?T S ^Z leSS than the m,ni7,urn of the ddta series or heater l^n 
the maximum of he oat, series. The curve corresponding to Equation (4) 

STuTto fJlSfS? *1 f f rel i abin ^ a Bol2ano "ethod t useo f r t 
fLitl w fV e E* Clear of the effect of an * ^flexions before the 

h?S rr pr^e a d Ph reliIbte r : t,0n ^ **" ""■ Th1s * ethod of «&2 

S» SlS&^pS2^ T "*" ' ike,ih00d est1 * ates * *• - fr ^ 

The Asymptotic Standard Error of Estimate 

KiSSiSS!*!! 1 varidnce of estimate equation the required partial 
derivatives are t putting z = lnx t p 



3z/3a = {t/3b^/6) - l/(9b 2 /3) ♦ fa*/*} 3 

»*/» - 3a{(t/3blA) . 1/(^/3) # b l/3} 2 { . t/(18b 7/6 ) «. 

2/(^7b 5 /3) ♦ l/(3b 2 /3)> 
3z/3m * 1 

and the inverse dispersion matrix V"- is 

and the inverse dispersion matrix V" 1 is 



D(9) 
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r i 



b/a 2 l/» 


1/a 2 


1/a tf (b ) 


l/(a(b-l)} 


1/a 2 l/{a(b-l)} 


a 2 /(b-2) 


the Triganuna Function 


ar*w 



lm I Z -„ () - n ? Vi 0) are eva1uated at the estimated values I and % 
with the appropriate t for the required return period. The elements of the 
matrix are evaluated at a and b and inverted to give the requ red var?ances 
and covanences Substitution in the general variance of est mate ewJtion 
Si r e V V ent ( ]? 1 " ^ ** As ^ totic Standard Error of i of ?he ? 



S,E. of x t = x t {var(lnx t )} 1 / 2 



! 
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Fitting by Moments 

l£\2ETlEi I.!* P,d - f ' g1ven by Equation UK re-arranging and 
replacing then, with their sample estimates gives 

* " (Se > ,/2 0(13) 

6 " ™^ 0(14) 

8 ■ X * 28/e > DOS) 
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where x, t and ^ --espectively are computed from Equations G(l), G(2) and 
G(3) with x replaced by its natural logarithm, lnx. 

Substitution of a\ ft and m^ in Equation 0(2) and using the appropriate value 
of t gives lnx t , and exponentiating gives the T ye&r event. 

Bobee, Reference (10) gives the variance of the T year event for a Pearson 
Type III distribution as a function of the second central moment, the skew 
coefficient and t. By a simple extension, for the Log-Pearson Type III, 



vaT 



where 

K w t + 



WM3 



and 



*1 6 63 l 63 Yl & T * " 6 6 "* 

Replacing >j 2 and Y t by their estimates, ^ and *g\ as previously computed 

for Equations 0(13) through D( 15) , and using the appropriate value of t, 

the variance of the logarithm of the T year event, var(lnx t ), is obtained. 

Then the Asymptotic Standard Error of Estimate of the T year event follows 
from 

S.C. of x t « x t {var(lnx t )} 1 / 2 
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D(18) 



D(19) 



-C17- 

.« 

REFERENCES 

(1) Nix, G.A., Condie, R., and Boone, L.G., Program HISTFLO, Flood 
Frequency Analysis Using Historic Information, Engineering Hydrology 
Section, Engineering Division, 1975. 

(2) Sangal, B.P., and Biswas, A.K., The 3 Parameter Log-Normal 
Distribution and its Applications in Hydrology, Water Resources 
Research, 6 (2) pp. 505-515, 1970. 

(3) Condie, R., and Kallio, R.W., Computer Programme ML3LN, Frequency 
Analysis, 3-Parameter Log-Normal Distribution, Parameters by Maximum 
Likelihood, Engineering Hydrology Section, Engineering Division, 
1975. 

(4) Bobee, B., The Log Pearson Type 3 Distribution and its Application 
in Hydrology, Water Resources Research, 11(5), 681-689, 1975. 

(5) Kimball, B.F., Sufficient Statistical Estimation Functions for the 
Parameters of the Distribution of Maximum Values, Annals of 
Mathematical Statistics, 17, 299, 1946. 

(6) Gumbel, E.J., Statistics of Extremes, p. 235, Columbia University 
Press, New York, 1960. 

(7) Kimball, B.F., An Approximation to the Sampling Variances of an 
Estimated Values of Given Frequency Based on Fit of Doubly 
Exponential Distribution of Maximum Values, Annals of Mathematical 
Statistics, 20, 110, 1949. 

(8) Kite, G.W., Frequency and Risk Analyses in Hydrology, Internal 
Document, Chapter 7.4.1, Water Resources Branch, Inland Waters 
Directorate, Environment Canada, 1976. 

(9) Wilson, E.B., and Hilferty, M.M., The Distribution of Chi-square, 
Proceedings of the National Academy of Science, U.S.A., 17, 684, 
1931. 

(10) Bobee, B., Sample Error of T-year Events Computed by Fitting a 
Pearson Type 3 Distribution, Water Resources Research, 9(5), 
1264-1270, 1973. 

(11) Programming Calcomp Pen Plotters, California Computers Products, 
Inc., 305 Muller Avenue, Anaheim, California, 92803, 1968. 



-D1- 



APPENDIX D 



COMPARISON OF FLOOD FREQUENCY USING LOG PEARSON TYPE III 
AND THREE-PARAMETER LOG NORMAL DISTRIBUTIONS 



Figure 0-1 GRAND RIVER AT SHAND DAM 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1978 
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Figure D-2 CONESTOGO KIVER AT CONESTOGO DAM 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 
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Figure 0-2 GRAND RIVER AT GALT 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 



Three-Parameter Log-Normal Distribution 
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Log Pearson Type III Distribution 
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Figure M GRAND RIVER AT BRANTFORD 
Flood Frequency Analysis 
Maximum Instantaneous Flow 194U- 1979 
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Figure D-i> SPEED RIVER BELOW GUELPH 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1950-1979 
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Figure D-6 NITH RIVtR AT NEW HAMBURG 

Flood Frequency Analysis 
Maximum Instantaneous Flow 1951-1979 
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Figure U- / NITH RIVER AT CANNING 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 
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APPENDIX E 

COMPARISON OF FLOOD FREQUENCY FOR THE GRAND RIVER AT GALT 1914-1979 AND 
1948-1979. 

The flood frequency analysis for the Grand River at Gait was carried out 
for the period 1948-1979 even though records extend back to 1914. 
Frequency relationships based on a consistent period of record at all 
stations were preferred for the development of frequency relationships at 
ungauged locations. 

Flood frequency analysis at Gait for the period of record 1914-1979 is 
presented in Tables E-l and E-2. 

In reviewing the flow records at Gait from 1914-1947, it was noted that: 

a) From 1914 to 1929 only maximum daily flows were available. Maximum 
instantaneous flows were estimated using the regression equations 
noted in Figure A-3. 

b) Data on reservoir operations at Shand Dam from 1942 to 1947 is 
sparse. The data is sufficient to provide rough estimates of the 
daily natural flows. The maximum instantaneous natural flows were 
estimated using the regression equations noted in Figure A-3. The 
greater error resulted in several inconsistencies between observed 
maximum instantaneous flows and the corresponding natural flows that 
were estimated. 

There is some evidence of an increasing trend in flood flows over the last 
66 years. The frequency analysis ignores this possible non-stationarity. 

A comparison of the frequency curves for 66 years of record (1914-1979)and 
the 32 years of record (1948-1979) is shown in Figure E-l. The comparison 
indicates that the shorter term analysis results in 12 to 14 percent higher 
flow estimates for the same return period than the longer term analysis. 
The shorter term frequency curve is within the upper 95% confidence 
interval of the longer term frequency curve. 
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TABLE E-2. NATURAL FLOti GRAND RIVER M SAL! 
MAXIMUM INSTANTANEOUS 1 LOW 1914 1979 
Jr! CUBIC FEET PER SEl ON© 

SAMPLI STATISTICS 
MEAN 2cH 23 •;, . D . 11376 . 8 

SAMPLE STATISTICS (LOGS) 
MEAN 10,0901 S.B. - , 41*12 

sAMPLfc KIM * 7210 SAMPLE hAX « 58700 
PARAMETERS LOR GUMBEL 1 A ■ .000111 
PARAMETERS FOR LOGNORMAL M 10.0901 S 
PARAMETERS FOR THREE PARAMETER LOGNORMAL A 

STA1 I S1ICS OF LOG<X -A) 
MEAN - 10.3216 S.D. - .3i»81 C.S. = .0335 C.K. - 

PARAMETERS FOR LOG PEARSON III BY MOMENTS A= .0633 B^ 48.5069 LOG(M>=13 . 1626 
PARAMETERS FOR LOG PEARSON 111 BY MAXIMUM LIKELIHOOD A ■ ".0710 B ■ 38.1182 
DISTRIBUTION STATISTICS- MEAN - 10.0901 S.D ■ .4383 C.S. * .3239 
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THREE PARAMETER 

LOGNORMAL 

FLOOD ST. ERROR 

ESTIMATE PERCENI 
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Figure E-l. GRAND RIVER AT GALT 

Flood Frequency Analysis 
of Natural Flows 



Maximum Instantaneous Flow (1914-1979) 

with 95% Confidence Limits 
Maximum Instantaneous Flow (1948-1979) 
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APPENDIX F 



EFFECT OF HURRICANE HAZEL ON THE FLOOD FREQUENCY CURVE 



The flood frequency analysis was carried out on maximum instantaneous flows 
generated during the spring. However, it is generally agreed that the 
upper end of the frequency curve will be dominated by events generated by 
tropical storms. Only one tropical storm has produced an annual maximum 
flow: Hurricane Hazel, October 1954. Figures D-l to D-5 show the effect of 
including this event on the frequency curves. 



Figure F-l. GRAND RIVEK AT bALT 

Flood Frequency Analysis 

Maximum Instantaneous Flow 1948-1979 



Hurricane Hazel not included. 



— — — Hurricane Hazel included. 
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Figure F-2. GRAND RIVER AT BRANTFORD 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1948-1979 



Hurricane Hazel not included. 



Hurricane Hazel included. 
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Figure F-3. SPEEU RIVER BELOW GUELPH 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1960-1979 



Hurricane Hazel not included. 



Hurricane Hazel included. 
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Figure F-4. NITH RIVER AT NEW HAMBURG 
Flood Frequency Analysis 
Maximum Instantaneous Flow 1951-1979 



Hurricane Hazel not included. 



— — — Hurricane Hazel included. 
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Figure F-S. NITH RIVER AT CANNING 

Flood Frequency Analysis 

Mdximuni Instantaneous Flow 1948-1979 



Hurricane Hazel not included. 



Hurricane Hazel included. 
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APPENDIX G 

COMPARISON OF FLOOD FREQUENCY WITH REGIONAL FLOOD FREQUENCY ANALYSIS 
COMPLETED BY ENVIRONMENT CANADA 

"Magnitude and Frequency of Flood in Southern Ontario", B.P. Sangal and 
R.W. Kallio, Technical Bulletin Series No. 99, Environment Canada, 1977 
details a regional flood frequency analysis for Southern Ontario. 

Southern Ontario was divided into nine hydrologic regions with roughly 
similar soils and physiographic characteristics. For each region a 
graphical relationship between the mean annual flood and drainage area was 
determined. Dimensionless frequency curves (ratio to mean annaul flood) 
based on the three-parameter log-normal distribution were derived for each 
region, giving the ratio of flood estimate probabilities. 

These curves for Region 4, Maitland-Thames-Grand System, have been 
reproduced as Figures E-l to 8. The corresponding results of this study 
are shown on the curves. 

In the case of mean annual flood versus drainage area, the higher flows 
resulting from the clay loam soils the north-western portion of the Grand 
River basin tend toward the maximum values of the regional curve; the lower 
flows from the gravel soils of the eastern Grand River basin tend toward 
the minimum values of the regional curve. 

In the case of the dimensionless frequency curves, the actual curve at all 
locations lies at or below the regional curve. 
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Figure G-l. Mean Annual Flood vs. 

Drainage Area for Region 4 
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Figure <l-2. Regional Frequency Curve for Region 4 
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Figure G-i- Regional Frequency Curve for Region 4 
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Figure b-4. Regional Frequency Curve for Region 4 
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Figure £-5. Regional Frequency Curve for Region 4 
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Figure G-6. Regional Frequency Curve for Region 4 
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Figure (j-7. Regional Frequency Curve for Region 4 
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Figure G-8. Regional Frequency Curve for Region 4 
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